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Abstract 

In this thesis, I look at the role of quantum entanglement in mixed-state quantum 
computation. The model we consider is the DQCl or 'power of one qubit' model. I 
show that there is minimal bipartite entanglement in a typical instance of the DQCl 
circuit and even put an upper bound on the possible amount of entanglement. The 
upper bound does not scale with the size of the system, making it hard to attribute 
the exponential speedup in the DQCl model to quantum entanglement. On the 
other hand, this limited amount of entanglement does not imply that the system 
is classically simulatable. This goes against the dictum that quantum evolutions 
involving states with limited amounts of entanglement can be simulated efficiently 
classically. A matrix product state (MPS) algorithm for a typical instance of the 
DQCl system requires exponential classical resources. This exposes a gap between 
the amount of entanglement and the amount of purely nonclassical correlations in a 
quantum system. 



X 

This gap, I suggest, can be filled by quantum discord. An entropic measure of 
quantumness that tries to capture the quantum notion of a system being disturbed 
by a measurement, it is a true measure of purely nonclassical correlations. I calculate 
it in a typical instance of the DQCl circuit and find that the amount of discord is 
a constant fraction of the maximum possible discord for a system of that size. This 
allows an interpretation of quantum discord as the resource that drives mixed-state 
quantum computation. I also study quantum discord as a quantity of independent 
interest. Its role in the phenomenon of entanglement distribution is studied through 
an easily comprehensible example. 

This thesis also contains discussions on the relation between the complexity 
classes P, BQP and DQCl. Additional material is presented on the connections 
between the DQCl model, Jones polynomials and statistical mechanics. The thesis 
concludes with a discussion of a few open problems related to the DQCl model, the 
quantum discord and their scope in quantum information science. 
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Chapter 1 



Introduction 



All great deeds, and all great thoughts have a ridiculous beginning. Great works 
are often born on a street corner or in a restaurant's revolving door. 

- Albert Camus 



The field of quantum information science is largely propelled by the discovery 
of quantum algorithms that promise unprecedented advantages over the best known 
classical algorithms. The algorithms of Shor, Grover and others are of more than 
just academic interest. They possess the potential of affecting our way of life in the 
information age. It is thus of economic, industrial, in addition to scientific interest 
to comprehend the physical principles these algorithms operate upon. These novel 
algorithms harness the intrinsic structure of quantum mechanics, believed to be the 
fundamental law of nature. Yet, by itself, this is not enough to completely unravel 
the mysteries of quantum computation. Part of this stems from the intrinsic mystery 
of quantum mechanics as a physical theory and its counterintuitive predictions when 
compared to daily life classical mechanics. For the remainder, there are at present 
certain gaps in our understanding of the working of quantum algorithms. It is some 
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of these gaps that we look at in this dissertation and attempt to plug. 

The ultimate aim of quantum information science is to build a functioning quan- 
tum computer. There have been spectacular advances in that direction. Multiple 
successful implementations of pure state algorithms in laboratories across the world 
have occurred. However, most of these schemes are limited by their lack of scalability. 
This is not entirely due to a lack of experimental capabilities. Quantum systems are 
intrinsically hard to isolate from their environment and manipulate. As the system 
gets bigger, the challenges mount. The higher the number of degrees of freedom, the 
higher the number of ways it can couple to the environment. Thus, if we are to build 
a functioning quantum computer, we must learn to operate with mixed states, the 
natural outcome of environmental decoherence. This has generically been tackled by 
quantum error correction, which aims at keeping quantum states as pure as possible 
through their quantum evolution. It was the notion of quantum error correction 
that first made the dream of a quantum computer realistic. Excellent techniques 
have since been devised to tackle an array of errors that might possibly plague a 
quantum computation, and this continues to be an active area of contemporary re- 
search. One needs to ensure that the error correcting part of the system is itself 
error-free, and most techniques of effective quantum error correction eventually rely 
on concatenation of several layers of error correction. The success of this scheme 
usually requisites the need to reduce errors below a threshold by successive layers 
of error correction, that is, fault tolerance for all operations. Though not entirely 
Sisyphean^, it is undeniable that the exercise of quantum computation aggregated 
with error correction makes the whole endeavor more challenging. 

A different approach would envisage a quantum computer using mixed states 
themselves. Their robustness against environmental decoherence and the reduction 
of some of the complications of quantum error correction is a worthy enough motiva- 

^Sisyphus, of Greek mythology, was punished by being cursed to roll a large boulder up 
a hill, only to watch it roll down again, and to repeat this throughout eternity. The same 
character lends his name to the Sisyphus method of laser cooling of atoms. 
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tion. Additionally, one might be interested in the computational potential of mixed 
quantum states by themselves. This is where the gaps in our theoretical understand- 
ing about quantum algorithms arrests our progress. Though some attempts have 
been made, most of which rely on distilling pure states out of them, e.g., magic state 
distillation [10], mixed-state quantum computation by itself is particularly poorly 
understood. In this thesis, our endeavor will be to address this situation. 

In 1998, E. Knill and R. Lafiamme presented the first intrinsically mixed-state 
scheme of quantum computation ["]. It was dubbed the 'power of one qubit' or 
DQCl (deterministic quantum computation with one pure qubit) model, as it re- 
quired only one pure qubit for its operation. The setup for this scheme is as follows: 
The top qubit is first acted upon by a Hadamard gate, given by 



|0)(0| -]H 
P 



l^= f^n 




Figure 1.1: The DQCl circuit 



H 



which effects the transformations 

The state of the system before the controlled unitary is 

1 



2„^„^, + |0)(l| + |l)(0| + |l)(l|)®p. 

Next the controlled gate means that the unitary f/„ acts on the lower set of n qubits 
in the state p if and only if the top (control) qubit is in the state |1). Otherwise it 
does not and /„ is applied. After this gate, the state of the system is 

^ = 1 (|0)(0| ^p + \0) (i| ® pui + |1) (0| ® U^p + |1)(1| ® f/„pf/t) . 
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This is the final state of the system and can compactly be written as 

^^\( p pul 

Suppose now that we make a measurement of the top qubit in the X basis, X being 
the Pauli operator 

.J- 

\l 

We leave the n lower qubits untouched. The expectation of this measurement is 
given, quite simply, by 



(X)= = tr[(X ® QE] = tr 



.Un + Ul 



Re[tr(pt/„)]. 



A measurement in the Y basis similarly gives an expectation value of Im[tr(pf/„)]. 
This is, in effect, the DQCl model of quantum computation. 

Suppose now that we prepared p in a pure state p = |$)($| that was an eigenstate 
of Un such that f/„ |$) = e*"^ |$). Then 

(X)= = cos^. 

Then L measurements of the top qubit in the X basis would provide us with an 
estimate of 0, via cos0 to within an accuracy of 1/vZ- This is, in effect, the phase 
estimation algorithm so commonly used in quantum algorithms. If, on the other 
hand, p = In/'^"', the completely mixed state, then L measurements will give us with 
an accuracy of 1/VL, the quantity 

{X)s = Re[tr(f/„)]/2^ 

This is the real part of a quantity called the normalized trace of the unitary matrix 
Un, and the scheme gives us an estimate for it with a fixed accuracy in a number 
of trials that does not scale with the size of the unitary matrix. Calculation of 
fidelity decay in quantum chaos, the evaluation of expectation values in condensed 
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matter physics and as shown very recently, even problems in quantum metrology, 
quadratically signed weight enumerators and estimation of Jones polynomials from 
knot theory and discrete quantum field theories can be reduced to the evaluation of 
normalized traces of particular unitary matrices. 

As unitaries conserve the purity of a state, tr(S^) = tr(p^). For p = In/'^"', 
the states involved in the DQCl model are thus very closed to being maximally 
mixed. Notwithstanding, this algorithm provides an exponential speedup over the 
best known classical algorithm. This shows that to harvest the benefits of exponen- 
tial speedups in quantum computation, one does not necessarily need to work with 
pure states. To be able to harness this quantum advantage optimally in the face of 
a decoherent environment, one needs to understand the very simple DQCl model 
better. Thus, the simple question I address in this thesis is the following- Why does 
the DQCl model work? This, I will attempt to answer in the subsequent chapters 
and as we will see, the answer is not simple. 

Quantum entanglement is the currency within the world of quantum information 
science. Over the last decade, it has come to be believed that entanglement is the 
resource that facilitates quantum advantages and speedups. This is based on evidence 
and knowledge that has been accrued since the early days of quantum information 
in the mid 1990's [48]. For instance, simulations have been done that study the rise 
and fall in the amount of entanglement in quantum evolutions that implement Shor's 
algorithm. The step having the maximum amount of entanglement is the so called 
'quantum Fourier transform', which has come be believed as the quantum soul of 
the whole algorithm. Quantum Fourier transform is a period finding technique, just 
like the conventional Fourier transform. The novelty, due to Shor, was to apply it 
over the field of integers, to find the period of the modular exponentiation function, 
to which integer factorization and discrete logarithm evaluation can be reduced ['2(i]. 
Grover's search algorithm has also been studied from an entanglement perspective. 
It has been found that the quadratic speedup in the algorithm is directly related to 
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the flow of entanglement in the algorithm. 

Preceding this line of work on quantum algorithms is the concept of quantum tele- 
portation [62]. Quantum teleportation is the ability to transfer the quantum state of 
a system from one point in space to another by just classical communication and a 
quantum resource. In its simplest incarnation, a state of the form |\1') = a |0) +/3 |1) 
can be teleported from Alice to Bob by sending just two bits of classical information 
if they already share a maximally entangled Bell state, say, 4=(|00) + |11)). Note 
that |\E') is completely denoted by three real parameters. This saving in the commu- 
nication cost stems from the inherent power of quantum mechanics. It is this power 
that the field of quantum information science strives to harness for problems more 
realistic than teleportation. Motivated by this, measures of entanglement are often 
defined in terms of the number of Bell states that are needed to produce or can be 
extracted from a certain state by local operations and classical communications. In 
2005, it was shown by Masanes that all entangled states (including bound entangled) 
can enhance the teleporting power of some quantum state [61]. 

Ever since the work on quantum teleportation, which is a genuinely quantum phe- 
nomenon, the research in quantum information science has been largely dominated 
by quantum entanglement. Evidence for this fact is the profuseness of papers begin- 
ning with the sentence - "Entanglement is a useful resource for quantum information 
and computation." Tragically, this realization does not advance much the cause of 
quantum information science. An instance inspiring such an opinion is the existence 
of quantum systems in highly entangled states that are computationally no more 
powerful that classical ones. Such a class of states are the so-called stabilizer states. 
Though highly entangled, their dynamics can be efficiently simulated classically. The 
twist in the story, is however, a result by Jozsa and Linden [~r2] which shows that 
for any pure-state quantum computation to be exponentially faster than its classical 
counterpart, the extent of entanglement in the system must not be bounded by the 
system size. 
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It is in the light of these circumstances that I analyze the DQCl model in Chapter 
3. The findings are intriguing. I found the first example of a unitary matrix for which 
the DQCl circuit showed a non-zero amount of entanglement. I then generalized 
this to create a family of unitaries which had the same behavior in terms of the 
entanglement content of the DQCl state. In collaboration, I then looked for the 
presence of bipartite entanglement in typical instances of the DQCl circuit. However, 
the amount of entanglement was small, in that the ratio of entanglement present to 
that maximally possible goes to zero in the limit of larger and larger systems. Since 
exponential speedups are determined asymptotically in the limit of large problem 
size, these results suggest that bipartite entanglement cannot be the reason behind 
the success of DQCl. I was also able to provide analytical upper bounds on the 
amount of entanglement (as measured by the negativity). This bound does not 
scale with size of the problem, and is the most crucial detriment in attributing the 
exponential speedup in the DQCl model to entanglement. In the course of my studies 
into the entanglement of the DQCl model, I was never able to find a unitary matrix 
which showed non-trivial negativity for a (polarization of the top qubit) less than 0.5. 
This presents a problem of the most singular nature. For a < 0.5 the DQCl state is 
at best bound entangled. The computational efficiency of the model, however, has 
no discontinuity at this value. If entanglement were behind the exponential speedup 
in the DQCl model, it would imply that this particular entanglement cannot be 
exploited for any other task (as it is bound and therefore cannot be distilled to form 
Bell states). This certainly would be the most peculiar state of affairs !! 

Another facet of the result by Jozsa and Linden is that quantum systems with 
small amounts of entanglement can be simulated classically. This is one of the reasons 
for believing entanglement to be a resource in the first place. This provided me with 
the impetus to study the classical simulat ability of the DQCl model. As presented 
in Chapter 4, I found that Matrix Product State (MPS) algorithms are incapable 
of simulating typical instances of the DQCl system efficiently in spite of its having 
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a limited amount of entanglement. The degree of correlations, as measured by the 
rank of the operator Schmidt decomposition, scaled exponentially with the size of 
the system. This I first showed by numerical simulation. Then, in collaboration with 
Guifre Vidal, I was able to provide an analytical proof that the operator Schmidt 
rank did indeed scale exponentially with the problem size. MPS techniques are the 
most commonly used for simulating quantum systems with limited amounts of entan- 
glement. My results, therefore, exclude a large family of algorithms from simulating 
the DQCl model. Simultaneously, it shows the claim made at the beginning of this 
paragraph to be invalid in the case of mixed quantum systems. 

The results of Chapters 3 and 4 leave us with a very uneasy realization. A typical 
instance of the DQCl model has very little entanglement, and yet it cannot be 
simulated classically. The natural conclusion is that there must be more to quantum 
computation that just entanglement. This is the scenario I investigate in Chapter 
5. I suggest that quantum discord^, a notion of quantum correlations more general 
than entanglement, is possibly that missing element which captures the power of 
a quantum computation completely. In particular, I show that across the most 
important split in the DQCl circuit, which is shown to be separable in Chapter 3, 
there is a non-zero amount of discord. The analytic calculations done by me were 
for the set of random Haar-distributed unitaries. This set is, however, not efficiently 
implementable in terms of any gate set universal for quantum computation, but 
sets of pseudo-random unitaries exist which are efficiently implementable and have 
eigenvalue statistics that mimic that of the Haar-distributed unitaries. I used these 
for numerical confirmation of my analytic results. Finally, the ratio of this discord 
to the maximum possible does not vanish in the limit of large systems. Thus I 
present quantum discord as the potential resource behind quantum advantages in 
computation. 

Quantum discord turns out to be an interesting and useful quantity in its own 



^This quantity [63, 44] was brought to my attention by Anil Shaji. 
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right, with apphcations wider than just to mixed state quantum computation. I have 
prehminary evidence that it plays a role in quantum communication, too, at least in 
the phenomenon of distributing entanglement. The details of this inference and the 
role of discord in entanglement distribution is presented in Section 2.3. 

In addition to studying the underpinnings of the DQCl model, I have studied 
an interesting connection between the DQCl model and knot theory. The DQCl 
model is equivalent to approximately evaluating the Jones polynomial of the trace 
closure of a braid. Rather, this problem is actually complete for the complexity 
class DQCl. The problem of approximately evaluating the Jones polynomial of 
the plat closure of a braid is BQP-complete. It is believed that the class DQCl 
actually lies in territory betwixt P and BQP, as shown in Fig 2.1. Since evaluation 
of the Jones polynomial is equivalent to evaluating the partition function of certain 
spin models like the Ising model, and different kind of braid closures correspond 
to different boundary conditions, this problem shows us a path for understanding 
the relations between complexity classes through physical models and theories like 
statistical mechanics. Thus, studying the DQCl model and its correspondence to 
Jones polynomials promises to be a very fruitful one, and I discuss this in a little 
more detail in Appendix A. 2. 

The above explorations by no means exhaust the class of problems a mixed-state 
quantum computer or the DQCl model can solve (See Fig 1.2 for some additional 
possibilities). One of the more exotic applications is through the quadratically signed 
weight enumerators (QSWEs). We will mention this as well in brief in Appendix A.l 
and not go into any details of it in this thesis. 
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Figure 1.2: A schematic that displays the rich interdisciphnary nature of quantum 
computation in general, and DQCl in particular. Black lines show the connections 
already known in the literature. Green lines are the connections first made in this 
thesis. The blue line show some of the possible connections that can be made and 
the potential for exploring the DQCl model from several independent points of view. 
'Trace' and 'Plat' refer to closures of braids whose Jones polynomials are of interest, 
while 'Periodic' and 'Aperiodic' refer to boundary conditions for Ising models. 'Graph 
enumeration' refers to the class of graph theoretic problems like the evaluation of 
the Tutte, the dichromatic, or the chromatic polynomial. The Kauffman bracket 
polynomial is a knot invariant which is a generalization of the Jones polynomial, and 
is related to partition functions of more esoteric statistical mechanical models like 
the vertex model and the IRF (interactions round a face) model. For further details, 
see Appendix A. 
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1.1 List of Publications 

It is unlikely that any physics PhD thesis in this age of instant electronic collaboration 
is a collation of entirely original results. So it is with this dissertation. In addition 
to the background material taken from the literature, several new results presented 
in this dissertation were arrived at in collaboration with other researchers. Several 
of these have been published in refereed journals. This dissertation contains only my 
research on mixed-state quantum computation and the DQCl model. In addition to 
this research, I have worked on two other topics during my graduate career at UNM, 
viz, quantum metrology and the abstract theory of entanglement. I provide below 
the list of my publications based on work done as a PhD student, classified according 
to which of the three topics they pertain to. For the papers on mixed-state quantum 
computation, I give which chapter each corresponds to in this dissertation. 

Mixed-State Quantum Computation 

1. Quantum discord and the power of one qubit, Animesh Datta, Anil Shaji, 
Carlton M. Caves, Phys. Rev. Lett, 100, 050502 (2008). (Chapter 5) 

2. Entanglement and the Power of One Qubit, Animesh Datta, Steven T. Flam- 
mia, Carlton M. Caves, Phys. Rev. A, 72, 042316 (2005). (Chapter 3) 

3. Role of entanglement and correlations in mixed-state quantum computation, 
Animesh Datta, Guifre Vidal, Phys. Rev. A, 75, 042310 (2007). (Chapter 4) 

Quantum Metrology 

4. Quantum Metrology: Dynamics vs Entanglement, Sergio Boixo, Animesh Datta, 
Steven T. Flammia, Matthew J. Davis, Anil Shaji Carlton M. Caves, Submitted 
to Physical Review Letters. 

5. Quantum-limited metrology with product states, Sergio Boixo, Animesh Datta, 
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Steven T. Flaminia, Anil Shaji, Emilio Bagan, Carlton M. Caves, Phys. Rev. 
A, 77, 012317 (2008). 

6. On Decoherence in Quantum Clock Synchronization, Sergio Boixo, Carlton M. 
Caves, Animesh Datta, Anil Shaji, Laser Physics: Special Issue on Quantum 
Information and Quantum Computation, 16, 1525 (2006). 

Quantum Entanglement 

7. Doubly constrained bounds on the entanglement of formation, Animesh Datta, 
Steven T. Flammia, Anil Shaji, Carlton M. Caves, arXiv:quant-ph/0608086. 

8. Constrained bounds on measures of entanglement, Animesh Datta, Steven T. 
Flammia, Anil Shaji, Carhon M. Caves, Phys. Rev. A, 75, 062117 (2007) 

The remainder of the Introduction explains in brief the results obtained in work 
not included in this dissertation. 

The study of entanglement in a mixed quantum state led to some work on the 
general theory of quantum entanglement, which though not a part of this dissertation 
is listed below. This project considered entanglement measures constructed from two 
positive, but not completely positive, maps on density operators that were used as 
constraints in placing bounds on the entanglement of formation, the tangle, and 
the concurrence of 4 x A^ mixed states. The norm-based entanglement measures 
constructed from these two maps, called negativity and phi-negativity, respectively, 
lead to two sets of bounds on the entanglement of formation, the tangle, and the 
concurrence. These bounds were compared, and we identified the sets of 4 x A^ density 
operators for which the bounds from one constraint are better than the bounds from 
the other. Extensions of our results to bipartite states of higher dimensions and with 
more than two constraints were also discussed. 

Another area of my research has been quantum metrology. The lure of beating the 
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shot-noise limit without using entanglement led my collaborators and myself to study 
various quantum clock synchronization protocols. We showed that the proposed 
schemes for beating the standard quantum limit without entanglement require as 
much resources as those using entanglement. We later turned to surpassing even 
the Heisenberg limit and proposed a scheme that involved product input states and 
measurements. Our results also provided evidence that the entanglement generated 
in the process was not responsible for the enhancement in precision. The next step 
involved proposing a demonstration of such a protocol using two-component Bose- 
Einstein condensates in ^^Rb. This is the first realistic laboratory proposal of beating 
the 1/n limit in metrology, certainly one using BECs 
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Chapter 2 



Background 



Take nothing on its looks; take everything on evidence. There's no better rule. 

-Charles Dickens in Great Expectations 



The intent of this chapter is to present research that is not entirely original, 
but nonetheless vital for the appreciation of the results appearing in the following 
chapters. Almost all of it can be found in the literature, and this will be pointed 
out at appropriate spots. This chapter will set the context and relevance of the 
DQCl model. Subsequently we will outline the basics of entanglement theory and its 
relevance in quantum computation. Later we will present arguments and motivations 
for considering a quantity such as quantum discord and its relevance to quantum 
information science. There is some original research in this chapter, in the later 
parts of Sec 2.2. All of Sec 2.3 is original, dealing with quantum discord in quantum 
communication. Though not a central theme of this thesis, this role is vital to our 
understanding of quantum discord. 

We begin this chapter with the requirements for a scalable quantum computation. 
We put the DQCl model in perspective with the rest of quantum computation. We 
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discuss the relevance of the DQCl model in several problems of physical importance 
and its place in computational complexity theory. We will briefly discuss the role of 
entanglement in all this and realize its shortcomings. Then we will move on to study 
a quantity that we think might be able to explain phenomena that entanglement 
fails at. This quantity, quantum discord, will be discussed in detail. We will explain 
how it is evaluated, and work out an example. Then we will study a phenomenon, 
called entanglement distribution, through a simple example and show, for the first 
time, the role of discord in its success. 

Without further ado, let us begin by what are postulated to be the requirements 
for building a fully controllable and scalable quantum computer. Known as DiVin- 
cenzo's criteria [22], they are as follows: 

1. Be scalable system in terms of well defined qubits (2-level quantum systems) 

2. Be initializable into a simple quantum state such as |00 ■ ■ ■ 000) 

3. Have decoherence times long compared to gate operation times 

4. Have a universal set of quantum gates 

5. Permit high efficiency, single qubit measurements 

However, fully controllable and scalable quantum computers are likely many years 
from realization. This motivates study and development of somewhat less ambitious 
quantum information processors, defined as devices that fail to satisfy one or more of 
DiVincenzo's five criteria for a quantum computer [22] listed above. An example of 
such a quantum information processor is a mixed-state quantum system, which fails 
to pass DiVincenzo's second requirement, that the system be prepared in an initial 
pure state. 

The prime example of a mixed-state quantum information processor is provided 
by liquid-state NMR experiments in quantum information processing [19]. Current 
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NMR experiments, which operate with initial states that are highly mixed thermal 
states, use a technique called pseudo-pure-state synthesis to process the initial ther- 
mal state and thereby to simulate pure-state quantum information processing. This 
technique suffers from an exponential loss of signal strength as the number of qubits 
per molecule increases and thus is not scalable. There is a different technique for 
processing the initial thermal state, called algorithmic cooling [ ], which pumps 
entropy from a subset of qubits into the remaining qubits, leaving the special subset 
in a pure state and the remaining qubits maximally mixed. Algorithmic cooling pro- 
vides an in-principle method for making liquid-state NMR — or any qubit system that 
begins in a thermal state — scalable, in essence by providing an efficient algorithmic 
method for cooling a subset of the initially thermal qubits to a pure state, thereby 
satisfying DiVincenzo's second criterion. 

Knill and Laflamme [55] proposed a related mixed-state computational model, 
which they called DQCl, in which there is just one initial pure qubit, along with 
n qubits in the maximally mixed state. Although provably less powerful than a 
pure-state quantum computer [-], DQCl can perform some computational tasks 
efficiently for which there are no known polynomial time classical algorithms. Knill 
and Laflamme referred to the power of this mixed-state computational model as 
the "power of one qubit." In particular, a DQCl quantum circuit can be used to 
evaluate, with fixed accuracy independent of n, the normalized trace, tr(f/„)/2", of 
any n-qubit unitary operator Un that can be implemented efficiently in terms of 
quantum gates [55, 57]. In Sec. 3.1 we consider whether there might be efficient 
classical algorithms for estimating the normalized trace, and we conclude that this 
is unlikely. 

The efficient quantum algorithm for estimating the normalized trace provides 
an exponential speedup over the best known classical algorithm for simulations of 
some quantum processes. It is known how to reduce the problem of measuring 
the average fidelity decay of a quantum map under small perturbations to that of 
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evaluating the normalized trace of a unitary matrix [ ]. The DQCl algorithm for 
this instance provides an exponential speedup. However, in testing the integrability 
of the classical limit of a quantum system, the same basic algorithm provides a 
square-root speedup [67]. Evaluating the local density of states, however, can be 
executed with an exponential advantage using only one pure qubit [i ]. There has 
also been an attempt to exploit the DQCl system in quantum metrology [8]. 

A very interesting application of the DQCl model has been to topological quan- 
tum computation. It has been shown that the evaluation of the plat closure of the 
Jones polynomial of the same braid is BQP-complete [31, 32, 33]. An explicit quan- 
tum algorithm for this task was presented in 2006 ['] in the language of quantum 
circuits. This technique has now been extended to more general topological invari- 
ants like the HOMFLYPT polynomial and the Kauffman polynomial [82]. Recently, 
it was shown that the DQCl model can be used to evaluate the Jones polynomial 
of the trace closure of a braid efficiently [51]. Indeed the problem of evaluating 
the Jones polynomial for the trace closure of a braid is complete for the complexity 
class DQCl. The notion of DQCl as a quantum complexity class was developed in 
[74]. Since then, one of the important questions has been the relation of this class 
to other well known classes like P and BQP. As has already been shown, DQCl is 
less powerful than BQP, a full scale pure state quantum computer [2]. Still, DQCl 
can efficiently simulate classical circuits of logarithmic depth [2], a class known as 
NCI, which is contained in P, the class of problems that can be solved on a classical 
polynomial in polynomial time. On the other hand, as we show in Section 3.1, it is 
unlikely that DQCl is contained in P. This leads to a very intriguing picture of the 
relative powers of these complexity classes in the polynomial hierarchy as shown in 
Fig 2.1. 

This figure has a number of intriguing question lurking underneath. The least of 
these is a proof of the location of DQCl in this figure. Also, if this figure is correct, 
as we believe at present, then DQCl is really a class that lies between quantum and 
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Figure 2.1: The relative hierarchy of the complexity classes P, BQP and DQCl. 



classical computation. One would also like to understand the physical principle that 
governs this relative positioning. Such a sort of understanding will tell us of the 
particular aspects of nature that help us beat classical computers using quantum 
mechanics. We have been pursuing this program with some success. We describe 
this is Appendix A. 

Study of the power of one qubit is motivated partly by NMR experiments, but 
our primary motivation in this thesis is to investigate the role of entanglement in 
quantum computation, using DQCl as a theoretical test bed for the investigation. 
For pure-state quantum computers, Jozsa and Linden [')'2] have shown that exponen- 
tial speedup over a classical computer requires that entanglement not be restricted to 
blocks of qubits of fixed size as problem size increases. Entanglement whose extent in- 
creases with problem size is thus a necessary prerequisite for the exponential speedup 
achieved by a pure-state quantum computer. On the other hand, the Gottesman- 
Knill theorem [62] demonstrates that global entanglement is far from sufficient for 
exponential speedup. While this means that the role of entanglement is not entirely 
understood for pure-state quantum computers, far less is known about the role of 
entanglement in mixed-state quantum computers. When applied to mixed-state com- 
putation, the Jozsa-Linden proof does not show that entanglement is a requirement 
for exponential speedup. Indeed, prior to the work reported in this dissertation, it 
has not previously been shown that there is any entanglement in the DQCl circuits 
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that provide an exponential speedup over classical algorithms. 

Evidently, to invoke the above results, one must first calculate the entanglement 
in a quantum computation. This exercise of calculating the entanglement in a general 
quantum state is highly non-trivial. Indeed, even the simpler task of proving the sep- 
arability of a general quantum state is NP-Hard [ ]. For a general mixed quantum 
state, the task is further hardened by the lack of a unique measure of entanglement. 
This is partly due to the fact that entanglement, instead of being a physical notion, 
is a mathematical construct. Though it has been identified with physical tasks that 
can be thought of as operational interpretations of entanglement, this has not sim- 
plified the task of actually calculating the numerical value of entanglement. The 
computational task is not eased by the definition of mixed-state measures through 
convex-roof extensions of pure state measures, most commonly the von-Neumann 
entropy of the reduced density matrix. One exception is the negativity, which is for- 
mally easy to define and computationally feasible to calculate, a rare feature for an 
entanglement measure. It is this measure that we use in this thesis, after elaborating 
upon it in the next section. For an extensive review on the theory of entanglement, 
the reader is referred to the recent article by the Horodecki family [IS]. 



2.1 Entanglement: Properties of negativity 

In this section we briefly review properties of negativity as an entanglement measure, 
focusing on those properties that we need in the subsequent analysis (for a thorough 
discussion of negativity, see Ref. [80]). 

Let A be an operator in the joint Hilbert space of two systems, system 1 of 
dimension di and system 2 of dimension d2- The partial transpose of A with respect 
to an orthonormal basis of system 2 is defined by taking the transpose of the matrix 
elements of A with respect to the system-2 indices. A partial transpose can also be 



Chapter 2. Background 20 

defined with respect to any basis of system 1. Partial transposition preserves tlie 
trace, and it commutes witli taking tlie adjoint. 

Tlie operator tliat results from partial transposition depends on which basis is 
used to define the transpose, but these different partial transposes are related by 
unitary transformations on the transposed system and thus have the same eigenvalues 
and singular values. Moreover, partial transposition on one of the systems is related 
to partial transposition on the other by an overall transposition, which also preserves 
eigenvalues and singular values. Despite the nonuniqueness of the partial transpose, 
we can talk meaningfully about its invariant properties, such as its eigenvalues and 
singular values. Similar considerations show that the eigenvalues and singular values 
are invariant under local unitary transformations. 



The singular values of an operator O are the eigenvalues of \J0'^0 = \0\ (or, 
equivalently, of vOO^). Any operator has a polar decomposition O = T\0\, where 
T is a unitary operator. Writing |0| = ly^S'iy, where W is the unitary that diago- 
nalizes \0\ and S is the diagonal matrix of singular values, we see that any operator 
can be written as O = VSW, where V = TW'^ and W are unitary operators. 

We denote a partial transpose of A generically by A. We write the eigenvalues 
of A as \j{A) and denote the singular values by Sj{A). If A is Hermitian, so is A, 
and the singular values of A, i.e., the eigenvalues of \A\, are the magnitudes of the 
eigenvalues, i.e, Sj{A) = |Aj(y4)|. 

If a joint density operator p of systems 1 and 2 is separable, its partial transpose p 
is a positive operator. This gives the Peres-Horodecki entanglement criterion [64, 45]: 
if p has a negative eigenvalue, then p is entangled (the converse is not generally 
true). The magnitude of the sum of the negative eigenvalues of the partial transpose, 
denoted by 

Afip)^- E ^^(Z^)' (2.1) 

A,(p)<0 
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is a measure of the amount of entanglement. Partial transposition preserves the 
trace, so tr(p) = 1, from which we get 

1 + 2A/-(p) = J2 \^M = tr|p| = M{p) , (2.2) 

j 

where M{p) is a closely related entanglement measure. The quantity A/"(p) was orig- 
inally called the negativity [ ]; we can distinguish the two measures by referring to 
Ai{p) as the multiplicative negativity, a name that emphasizes one of its key prop- 
erties and advantages over A/'(p). In this thesis, however, we use the multiplicative 
negativity exclusively and so refer to it simply as "the negativity" . 

The negativity Ai{p) equals one for separable states, and it is an entanglement 
monotone [bij], meaning that (i) it is a convex function of density operators and 
(ii) it does not increase under local operations and classical communication. The 
negativity has the property of being multiplicative in the sense that the Ai value 
for a state that is a product of states for many pairs of systems is the product of 
the Ad values for each of the pairs. By the same token, \ogAi{p), called the log- 
negativity, is additive, but the logarithm destroys convexity so the log-negativity is 
not an entanglement monotone [^'^]. For another point of view on the monotonicity 
of the log- negativity, see Ref . [(i5] . 

The minimum value of the negativity is one, but we need to know the maximum 
value to calibrate our results. Convexity guarantees that the maximum value is 
attained on pure states. We can find the maximum ["' ] by considering the Schmidt 
decomposition of a joint pure state of systems 1 and 2, 

d 

l^) = Ev^l^'^')' (2.3) 

where d = min((ii, ^2). Taking the partial transpose of p relative to the Schmidt 
basis of system 2 gives 

d 

P=Y1 V^^J^^k\j,k){k,j\ , (2.4) 

j,k=i 
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with eigenvectors and eigenvalues 

|j,j) , eigenvalue /ij, 
-7^i\j,k) ±\k,j)) , eigenvalue ±^/I^, j < k. (2.5) 



This gives a negativity 



d , d X 2 

j<k j,k=l ^1=1 ^ 



(2.6) 



The concavity of the square root implies ^ ■ ^/JIj < vd, with equality if and only if 
fij = 1/d for all j, i.e., \tp) is maximally entangled. We end up with 

l<Mip)<d. (2.7) 



2.2 Beyond Entanglement: Quantum Discord 

Characterizing and quantifying the information-processing capabilities offered by 
quantum phenomena like entanglement, superposition, and interference is one of the 
primary objectives of quantum information theory. Amongst these, entanglement has 
defined much of the research in quantum information science for almost a decade now. 
In spite of some progress [27, 52, 53], the precise role of entanglement in quantum 
information processing remains an open question [6, 9, 77, 54, 18]. It is quite well 
established that entanglement is essential for certain kinds of quantum-information 
tasks like teleportation and super-dense coding. In these cases, it is also known that 
the quantum enhancement must come from (pre-shared) entanglement between parts 
of the system. It is not known, however, if all information-processing tasks that can 
be done more efficiently with a quantum system than with a comparable classical 
system require entanglement as a resource. Indeed, there are several instances where 
we see a quantum advantage in the absence or near absence of entanglement. One 
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is quantum cryptography [5], where many protocols involve quantum states that are 
not entangled. These quantum cryptography protocols are provably more secure and 
thus better than the best known classical cryptographic techniques. This certainly 
cannot be attributed to entanglement. The second example we present is that of the 
DQCl model, the theme of this dissertation. As this thesis shows, in Chapter 3, the 
system has very little entanglement - and even that little vanishes asymptotically - 
yet it provides an exponential speedup. These examples drive home one point: all of 
quantum information science cannot be reduced to a study of quantum entanglement. 

This however leaves us in a lurch. What then explains all the quantum advan- 
tages ? Certainly, it has something to do with the structure and inherent non-locality 
of quantum mechanics. Interestingly, it is known that non-locality and entanglement 
are not equivalent features [j i]. Entanglement is the feature that having complete 
information about parts of a quantum system does not imply complete information 
of the whole system. This characterization goes all the way back to Schrodinger 
[69]. Not surprisingly, this is not the only characterization of a quantum system. 
There are other ones. For instance, the collapse of one part of a subsystem after 
a measurement of another is another feature unique to quantum systems. One can 
conceive of defining a quantity that captures this feature. Such a quantity is the 
quantum discord [63]. 

Quantum discord has been related to the superposition principle and the vanish- 
ing of discord shown to be a criterion for the preferred, effectively classical states of 
a system, i.e., the pointer states [63]. In addition, it has been used in analyzing the 
powers of a quantum Maxwell's demon [84] and also in the study of pure quantum 
states as a resource [16]. 

We start with a discussion of quantum discord, its definition and its relevance in 
quantum information theory. Consider the following two-qubit separable state 



1 
P=4 



(l+)(+l®|0)(0|) + (|-)(-|®|l)(l|) + (|0)(0|®|-X-|) + (|l)(l|®|+)(+|) , (2.8) 
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in which four nonorthogonal states of the first qubit are correlated with four nonorthog- 
onal states of the second qubit. Such correlations cannot exist in any classical state 
of two bits. The extra correlations the quantum state can contain compared to an 
equivalent classical system with two bits could reasonably be called quantum correla- 
tions. Entanglement is certainly a kind of quantum correlation, but it is not the only 
kind. In other words, separable quantum states can have correlations that cannot 
be captured by a probability distribution defined over the states of an equivalent 
classical system. 

Quantum discord attempts to quantify all quantum correlations including entan- 
glement. It must be emphasized here that the discord supplements the measures 
of entanglement that can be defined on the system of interest. It aims to capture 
all the nonclassical correlations present in a system, those that can be identified as 
entanglement and then some more. 

The information-theoretic measure of correlations between two systems A and B 
is the mutual information, I{A : B) = H{A) + H{B) - H{A,B). If A and B are 
classical systems whose state is described by a probability distribution p{A, B), then 
H{-) denotes the Shannon entropy, H{p) = — YlijVj logPii where p is a probability 
vector. If A and B are quantum systems described by a combined density matrix pab-, 
then H{-) stands for the corresponding von Neumann entropy, H{p) = — tr(plogp). 

For classical probability distributions, Bayes's rule leads to an equivalent expres- 
sion for the mutual information, I{A : B) = H{B) — H{B\A), where the conditional 
entropy H{B\A) is an average of Shannon entropies for B, conditioned on the alter- 
natives for A. For quantum systems, we can regard this form for T{A : B) as defining 
a conditional entropy, but it is not an average of von Neumann entropies and is not 
necessarily nonnegative [ ]. 

Another way of generalizing the classical conditional entropy to the quantum case 
is to recognize that classically H{B\A) quantifies the ignorance about the system B 
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that remains if we make measurements to determine A. When A is a quantum 
system, the amount of information we can extract about it depends on the choice of 
measurement. If we restrict to projective measurements (a priori, we could consider 
POVMs, but we will show later that rank one projectors suffice) described by a 
complete set of orthogonal projectors, {Hj}, corresponding to outcomes j, then the 
state of B after a measurement is given by 

PB\j = tTM{lljPAB^j)/Pj, Pj = tTA,B{pAB'n.j). (2.9) 

A quantum analogue of the conditional entropy can then be defined as H{ji.y{B\A) = 
J2jPjHiPB\j) > 0. Since pb = J2jPjPb\j, the concavity of von Neumann entropy 
implies that H{B) > H^u }{B\A). We can now define an alternative quantum version 
of the mutual information, 

J{u,}{A : B) = H{B) - H{u,}{B\A) > 0. (2.10) 

Performing projective measurements onto a complete set of orthogonal states of 
A effectively removes all nonclassical correlations between B and A. In the post- 
measurement state, mutually orthogonal states of A are correlated with at most as 
many states of -B. It is easy to see that these sorts of correlations can be present in 
an equivalent classical system. If on the other hand, the system B is bigger than A, 
the correlations are certainly non-classical, or quantum. 

The value of J7{n }(^ : B) in Eq. (2.10) depends on the choice of {IIj}. We 
want J'{u }{A : B) to quantify all the classical correlations in pab, so we maximize 
J'^Yi}{A : B) over all {IIj} and define a measurement-independent mutual informa- 
tion J{A : B) = H{B) - H{B\A) > 0, where H{B\A) = mm^u,} Y.jPjH{pb\j) is a 
measurement-independent conditional information. Henderson and Vedral [44] inves- 
tigated how J{A : B) quantifies classical correlations. The criterion they postulated 
for a measure of classical correlations C were [44] 

1. C = for product states 
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2. C is invariant under local unitary transformations. This is because any change 
of basis should not affect the correlation between two subsystems. 

3. C is non-increasing under local operations. If the two subsystems evolve inde- 
pendently then the correlation between them cannot increase. 

4. C = H{pa) = H{pb) for pure states. 

One can then define a measure of purely quantum correlations as the difference of 
the total correlations in a system and C. Taking the quantum mutual information to 
be a measure of total correlations in a system, the purely quantum correlations can 
be measured by 

V{A, B) = I{A : B) - J{A : B) = H{B\A) - H{B\A). (2.11) 

This quantity was first called quantum discord by Ollivier and Zurek in [63] . Another 
operational approach to identifying the nature of correlations is presented in [37] . The 
discord of the state (2.8) is | log | = 0.311. 

This thesis is the first to propound a role for quantum discord in quantum com- 
putation and information. It is, therefore, imperative that we present all the known 
properties of discord here. This is the content of the next section. We will present 
several properties of quantum discord, some of which are present in prior literature. 
We will provide appropriate indications at requisite spots. 

2.2.1 Properties of Quantum Discord 

The first property that we present is that the discord is nonnegative and is zero 
for states with only classical correlations [14, 63]. Thus a nonzero value of V{A,B) 
indicates the presence of nonclassical correlations [63] . 

Theorem 1. Quantum discord is always positive, i.e., V{A,B) > 0. 
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Proof: Consider the joint state pab subject to one dimensional orthogonal mea- 
surements Uj = \ej){ej\ on B, extended to arbitrary (at most dim{BY) dimensions. 
Then 

Pj PA\J = tTBipAB'n^j) = {ej\pAB \ej) , Pj = tlBips'^^j) = {^j] pAB \ej) ■ 

Note that the measurement is made on the system B, while in the definition of 
discord, it was on A. Discord is not symmetric under the exchange of the subsystems, 
but this is not a concern as we just as well have proved the result for 'D^B, A). 

Suppose now that a system C interacts with B so as to make the desired mea- 
surement (f/ \ej) (g) |0) = \ej) (g) \fj)), leaving the state 

p'abc = Yl (^j\P^B \ek) ® \ej) {ek\ ® \fj) {fk\ ■ (2.12) 

If the eigendecomposition of pab = J2i '^i\^i){''"i\^ then 

p'abc = Y -^/(Ia, ej\ri){ri\lA, e^) (g) \ej) {ek\ ® \fj) {fk\ 

i,k,l 

I 

whereby 

Hip'ABc) = H{pab). 
Also, from Eq. 2.12, 

p'ab = 'YPjPAij^lejXcjl, so H{p'^b) = H{p) + YpjH{A\j), (2.13a) 
j j 

p'bc = Y\(i3){(ij\pB\ek){ek\®\fj){fk\, so H{psc) = H{pb), (2.13b) 

p'b = J^PjIeiXcjl, so H{pb) = H{p). (2.13c) 

j 

Now use the strong subadditivity of the von- Neumann entropy which is 

H{p'abc) + H{p'b) < HipUs) + Hip'^c)- (2-14) 
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Eqs. 2.13 reduces this to 

H(pab) + H(p) < H(p) + J2pjH(A\j) + H{pb), (2.15) 

j 
whereby 

H{n,M\B) ^ Y.P^HiA\j) > H{pab) - H{pb) = H{A\B). (2.16) 

3 

This, being true for all measurements, also holds for the minimum. So 



V{A,B) = min H{u,}iA\B) - H{A\B) > 0. 



D 



Having proved that the quantum discord is always nonnegative, it is worthwhile 
to seek the condition for it being zero. The reason is that the set of states with zero 
discord is exactly those which have no nonclassical correlations in it. This is the 
aim of the next theorem. It is evident that the condition for zero discord can be 
reduced to that of the equality in strong subadditivity in Eq 2.14. To that end, we 
will employ a result of Hayden et al. [ _] which we present below for completeness. 

Lemma 1 ([ ])• A state p'^Bc on7iA®'HB®'Hc satisfies strong subadditivity (Eq. 
2. 14) with equality if and only if there is a decomposition of system B as 

j 
into a direct sum of tensor products such that 

p'abc = ® QjPAB^ ® PABf 

j 

with states Pab^ on 1-La®'Hbl and Pab^ on TioR^'Hc, and a probability distribution 

i i i j 

Theorem 2. V{A, B) = if and only if the state pab is block diagonal in its own 
eigenbasis, that is 

Pab = 2_^ PjPabPj 
j 
where pab = J2j '^jPj' '"^^^^ {''"} ^ probability distribution. 
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Proof: The decomposition of the Hilbert space of B can be written as 

a a 

and nall/g = Sa/3^a- In our case, the state p'abc i^ invariant under the exchange of 
B and C relative to the measurement basis, here denoted by \Eaj) and \Faj). Thus 
using Lemma 1, we conclude that it must have the form 

P'aBC = 2^ QaPA\a ® PbC^ 
a 

where p%Q = IlaP%cIla, with Ua being projectors of the form 

Hq, = 2_^ \Eaj){Eaj\ ® \Eaj){Faj\. 



Thus 



and 



PbC = Z^P'jk \Eaj) {Eak\ ® \Faj) {Fak 
j,k 



P'aB = Yl ^°^^l" ®Pb = Yp'jj^lc'PAla ® \Eaj) {E^jl ■ 
a a,j 

Undoing the measurement, U \ej) |0) = \ej) \fj), gives 

PAB = (Ocl U^p'abcU |0c) = XI ^^PAW ® Pb- 

a 

Diagonalizing p% = ^j A"|A")(A°|, we get 

pAB = Y>^"qapA\a\>^'j){>^'j\- 

Relabelling, we have that the discord is zero if and only if 

PAB = YpjPA\j(^\Xj){Xj\ (2.17) 

j 

in the basis that diagonalizes pb- The a subspaces take into account that if the 

states pA\j are the same for different j, then we can attain zero discord by using any 

measurement in the subspace spanned by those values of j. 
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Diagonalising pA\j = J2kf^jk\f^jk){fi'jk\, we get that a state has zero discord if and 
only if 

Pab = ^ PjfXjklPjk, \){Pjk, \ I • (2. 18) 

jk 

Thus the eigenbasis of pab has a tree product structure Ipjk) ® |Aj). From Eq. 2.18, 
it is also evident that a state has zero discord if and only if it is block diagonal in its 
eigenbasis, that is, 

V{A,B) = {] iff pab = Y.P,PabPj (2.19) 

i 

with 

Pj = 2_^ IPjk, ^j)\pjk, Aj|. 
k 

which is the statement of the theorem. D 

It is now a valid question to ask for the maximum possible value of discord. The 
next theorem addresses just this point. 

Theorem 3. The value of quantum discord is upper bounded by the von-Neumann 
entropy of the measured subsystem, i.e., V{A,B) < H{A). 

Proof: Start with the eigendecomposition 

Pab = y^PgHg, 

a 

which yields 



ij' 



PjPB\j = ^PaPj\apB\a,: 
a 

where 

PB\a,j = trA(n.jPB'n-a)/Pj\a, 

is a pure state of 5. It follows from the pure-state decomposition 

PB\j = ^Pa\jPB\a,: 



a 
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that H{{a}\j) > H{pB\j). Thus 

HiAB) = Hi{a})>Hi{a}\j) 

= J2P^H{{a}\j) 
j 

> H{n,}{B\A) 

> H{B\A), (2.20) 

from which the marginal entropy H{A) follows as the upper bound on discord. D 

When the joint state pab is pure, H{A, B) and H{B\A) are zero, H{B) = H{A) = 
—H{B\A), and the discord is equal to its maximum value, H{A), which is a measure 
of entanglement for bipartite pure states. In other words, for pure states all nonclas- 
sical correlations characterized by quantum discord can be identified as entanglement 
as measured by the marginal entropy. 

Coming back to the question of mixed quantum states, we find that the actual 
evaluation of quantum discord involves the minimization of an entropic quantity over 
the set of all POVMS, as shown in Eq 2.11. This minimization, though perhaps easier 
than the one involved in the convex roof optimization needed for certain measures of 
entanglement, is nevertheless not trivial. Some simplification is in order as we show 
in the following two theorems. 

Theorem 4. The quantity V{A,B) = H{B\A) — H£{B\A) is a concave function on 
the setofPOVMsS. 

Proof: We start by stating that for quantum states the von Neumann entropy is 
concave, that is, 

H{Xp, + (1 - A)p2) > XH{pi) + (1 - X)H{p2) (2.21) 

for all < A < 1. To study the concavity of the classical measure H{B\A), let intro- 
duce Positive Operator Valued Measurements (POVMs) given by a set of quantum 
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operations 

A- = E^^°®^1' (2.22) 

a 

with 

J2aIA,^ = E, and 5^i?,=I 

a j 

Not confusing the POVM elements Aja with the system A on which they are exe- 
cuted, the post-measurement state of B is 

PB\J = tr^(^j(pAB))/pi = 5]]trA(AjaPAB^]«)/pi 

a 

= tr^ {pAB Y^ A]^Aja ) /Pj 

a 
= tTA{pABEj)/pj. (2.23) 

The post-measurement state is thus completely determined by the POVM. Now let 
us call the POVM S = {Ej} and define the functional F{pB,S) as 

FipB,S) = HiB\A) = Y,pMpb\,). (2.24) 

3 

Let C and T> be two POVMs and £ = \C + {1 — \)V their combination. S is certainly 
a POVM since the set of POVMs is convex. This is, in effect, tantamount to 

E,=\C, + il-\)D,. 

Then it is simple to show that 

Pej Ps\e, = A pcj pB\Cj + (1 - A) Pd, Pb\d, , Pe, = A pcj + (1 - A) po^ ■ 

Using these two relations and the fact that the von-Neumann entropy is concave, it 
can be seen that 

F{pB,S) = Ype^H{pb\e,) 

j 

> Yl ^Pc,H{pb\c^) + (1 - ^)Pd,H{pb\d,) 

j 
= \F{pB,C) + {l-\)F{pB,V). (2.25) 
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This shows that the objective function which is to be minimized to obtain the discord 
is concave over the convex set of POVMs . D 

Thus the minima that provides the value of discord will always be attained at the 
extremal points of the set of POVMs. The extreme points of the set of POVMs on 
a D dimensional system is no more than D^ dimensional; a necessary and sufficient 
condition for its extremality being that the eigenvectors of the POVM elements 
be linearly independent as operators [ ]. It consists of projectors of all ranks. 
Fortunately, the next theorem will provide a considerable simplification. 

Theorem 5. The minimum of the quantity T'(v4, B) = H{B\A) —H£{B\A) is always 
attained for a rank-1 POVM. 

Proof: We start by supposing a POVM on system A, whose elements can be fine- 
grained as 

Ej = ^Ejk- 

k 

Then 

PjkPBlJk = tTAipAsEjk), Pjk = tT{pABEjk). 

Evidently, YlikPjk = Vj whereby we can define pk\j = Pjk/Pj- Also, 

pB\j = ti a{pabEj)/pj = ^ -^tr A{pABEjk)/pjk = ^Pk\jpB\jk- (2.26) 

k ^^ k 

Now, 



J2PjH{pB\j) = J^Pi^ [J2p''\^P^\^^J 

j j \ k / 

> ^P3Pk\jH{pB\jk) 



= Y.p,kH{pB\jk). (2.27) 

Since any POVM element can be written in terms of its eigendecomposition, the 
minimum conditional entropy, and therefore the discord is always attained on a 
rank-1 POVM. D 
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A rank-1 POVM is extremal if and only if its elements are linearly independent 
[16]. Even after these simplifications, the minimization over all rank one POVMs 
remains a hard task for general Hilbert spaces. We are not aware of the exact 
computational complexity of the problem, which might be formidable, but one can 
envisage casting this as a semi-definite program. This particular aspect is still under 
investigation by us. Accordingly, this thesis deals only with instances where the 
minimization is done over a two dimensional Hilbert space. 

In light of the present scenario, we conclude the discussion on the evaluation of 
quantum discord by explicitly working out an example. The state we consider is 
the well known 2x4 bound entangled state due to Horodecki [li], written in the 
standard product basis as 
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for < p < 1. (2.28) 



/ 



That this state is bound entangled for all < p < lean be concluded from its 
positive partial transpose, and its entanglement is hidden in a very subtle way. The 
calculation of the discord begins with the evaluation of the von- Neumann entropy of 
p. Its spectrum is 



A 



P 



2p 2p 



°'°'°'l + 7p'l + 7p'l + 7p' 
l + 9p+ Up^ 



l + 12p + 23p2 - 70p3 + 98p4 



2(1 + Up + 4V 



l + 9p+ Up'^ 



1 + I2p + 23p^ - 70p'^ + 98p^ 



2(1 + 14p + 49p2) 



(2.29) 
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from which H{p) = H{X) = — tr(Alog(A)). The reduced density matrix of the 
smaller subsystem is 



J^ 



(2.30) 



l+7p 



whose von-Neumann entropy is trivially obtained. As is clear by now, we have 
chosen to make our measurement on the 2 dimensional subsystem. These projectors 
are generally defined as 



Hi 



Jo + a ■ (7 



h 



n. 



lo — a ■ (7 



I4, 



2 ' 2 

with Oi = sin ^ cos 0,02 = sin ^ sin 0,03 = cos^. Using these, we calculate 



(2.31) 



Pi = tr(nipni) 
P2 = tr(n2pn2) 



1 + 7p — cos 6 + p cos 6 
2 + 14j9 ' 

1 + 7p + cos 6 — p cos 6 
2 + 14p ' 



(2.32a) 
(2.32b) 



and 



Pi 



trM(nipni 



1 



X 



(2.33) 



Pi 4(1 + 7p)pi 

f l + 3p-cose{l-p) 2psmee'^ 

2p sin 6'e~*'^ 4p 2p sin 6'e*'^ 

2p sin 9e-''^ Ap 

\ x/l-p2(l-cos^) 2psin^e-^'^ 1 + 3p - cos^(l - p) / 



l-p2(l_coS^) \ 



2p sin Oe'^'^ 



1 



P2 



trM(n2pn2) 

P2 4(1 + 7p)p2 

I 1 + 3p + COS 6'(1 - p) -2p sin 6'e^'^ 
— 2psin6'e~*'^ Ap 

-2psin^e-*'^ 





(2.34) 



l-p2(i + cos^) \ 

— 2psin6'e*'^ 



\ a/1 -p2(l + cos( 





— 2psin^e*"^ 

Ap 
-2psin^e-*'^ 1 + 3p + cos^(l - p) / 
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Figure 2.2: The green dashed hne is the value of the discord when the projector 
(Eq 2.31) is parameterized by 9 = 7r/2, = 0. The blue dot-dashed line is when the 
projector is given by 6' = 0, = 0. The solid red line, which is the smaller of the 
two, is the quantum discord in the 2x4 bound entangled Horodecki state of Eq 2.28 
obtained by explicit numerical minimization. The expression for the discord is given 
by Eq 2.35. We see how the optimal projector changes from the former to the latter 
at p = 1/7. This exemplifies how the optimal measurement can change even for the 
same family of states and therefore is tricky to obtain analytically. 



This state provides a very good example when optimizing over the POVM though 
seemingly tractable analytically, actually leads to erroneous conclusions. This is due 
to the fact that the optimal measurement that minimizes the conditional entropy 
changes with the value of p. We can however do it numerically without much ado. 
The value of discord in this case is 

Ap 



V = Ho 



+ tr(A log(A)) + mill [piif (pi) + P2Hip2)] 

i + ipj (6»,0) 



(2.35) 



7^2 (■) being the binary Shannon entropy. We plot the result of the numerical mini- 
mization as the green line in Fig 2.2. The plot was obtained in a couple of minutes 
on a 1.83 GHz Dell XPS laptop. In addition to demonstrating the feasibility of the 
minimization in small instances. Fig 2.2 shows one very important feature. This 
state is bound entangled, meaning it has entanglement that cannot be distilled [ 16] . 
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Quantifying the amount of entanglement in a quantum state is generally hard, harder 
still for bound entangled states. As we have now shown, quantum discord can be 
used to used to quantify non-classical correlations in such states, which further makes 
discord a useful tool in the study of quantum states. Additionally, we are going to 
demonstrate an analytical technique that might be useful in the minimization as well 
its caveats. Let us now proceed towards that. 

It is worthwhile to point out at the outset that the optimal measurements that 
give us the actual value of discord are not given by the eigenvectors of the reduced 
density matrices. Indeed, an independent measure of non-classicality has been pro- 
posed where measurements are made on both the subsystems, and that too in the 
eigenbasis of their respective subsystems [ ]. In the minimization of quantum dis- 
cord, the hard part is often using the explicit analytical expressions of the eigenvalues 
of the states in Eqs 2.33, 2.34. Here is a handy trick that is often helpful. Consider 
p = 1. In that case the Eqs 2.33, 2.34 simplify drastically to 
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(2.36) 



and 
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(2.37) 
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both of which have spectrum 



A(pi 



A(pi 



4- v^e- 2^5 sine 4 + a/6 - 2y5sin0 



16 ' 16 

4- \/6 + 2^5 sine 4 + \/6 + 2A/5sine 



16 



16 



}■ 



(2.38) 



Simuhaneously, Pi = P2 = 1/2. Then the conditional quantum entropy is just 
H{X{pi)), the Shannon entropy of the vector in Eq 2.38, which is minimized when 
the distribution is most disparate, that is, 9 = 7r/2. We note that 0, having dropped 
out of the problem, can be set to zero. We now have at least one point where we 
have carried out the minimization explicitly. One might now be tempted to believe 
that this point gives the minimum for all p, and we would have succeeded in the 
minimization without doing it explicitly. Fig 2.2 shows the catch. The red line is the 
expression for discord obtained when setting 6 = 7r/2, = 0. We see that it matches 
with the actual minima for a range of p, but not all of it. The region where they 
diverge is p G {0, 1/7}. Interestingly enough, in this region, the optimal choice of 
the projector is given by ^ = 0, = 0. This is shown by the blue line, which we see 
lies exactly on top of the green line obtained by explicit numerical minimization in 

pe {0,1/7}. 

To reemphasize the point made at the beginning of the last paragraph, we present 
the eigenvectors of E(pi) and E(p2) of pi and p2 in Eqs. 2.36 and 2.37 respectively. 
For all values of 6 and 0, 
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V5+1 i(l> 
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V5+1 i(t> 
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(2.39) 
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and 



E(P2 



/ ^^Siip \ 
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/5-l^id 
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V5+1 2i<j, 
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\] 
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(2.40) 



and the ordering of the eigenvectors correspond to the eigenvalues in Eq 2.38. As is 
evident, these is no 6 dependence in these vectors. Accordingly, it is no surprise that 
they cannot provide any information about the optimal measurement that appears 
in the calculation of discord. 

As mentioned before, and reiterated by this example, there is a gap between 
discord and entanglement. Certainly, zero discord states are separable. The converse 
is, in general not true for mixed quantum states. In light of this, a valid question is 
whether there is a bound on the value of discord for separable states. This result can 
be used indirectly to detect an entangled state. In the following section we inquire 
about the value of quantum discord for separable states. 

2.2.2 Discord for separable states 



Let us consider a separable state 



PAB = ^Ptpf®P\ 



B 

i 1 



whereby 



B 



PA = ^Pipt and pB = ^Pip. 

i i 

The quantum discord of a quantum state is defined as 



V = H{pa) - H{pab) + min^ g^ H{pb\j), 

{n,} . 
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where 

tl Ai'^jpAB'^j) ^ / TT \ 

PB\j = ^ -, qj = tTA,B[PABi^j)- 

'ij 

It is known that for a separable state 

H{pab)>H{pa) and H{pab) > H{pb), (2.41) 

whereby 

V = H{pa) - H{pab) + mm ^ q^ H{pB\j) 

in, I ^ 

- '%^'^^^3H{pB\j) 

in,} . 
^ '^f^H\y2q,PB\y (2.42) 

where the last inequality follows from the concavity of H. Now 

Y^ QjPBlJ = Yl ^^^ O^jPAB^lj) = Y^ tTA (Jlj Pab) 
3 J J 

= tTA({J2^j)pAB]=pB. (2.43) 

J 

Therefore, we have, 

V < H{pb). 

Of course, one can choose not to invoke the separability condition first, and rather 
say that (since JT" > 0) 

V < H{pa) + H{pb) - H{pab), (2.44) 
and then employ (2.41) to conclude 

V<mm{H{A),H{B),I{A: B)), (2.45) 

which is probably a more sensible upper bound for separable states. This, in fact, 
leads to the conclusion that X> = for pure separable states. 
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2.3 Quantum discord in quantum communication 



The centerpiece of this thesis is undoubtedly mixed-state quantum information. 
Nonetheless, we have an instance from quantum communication where quantum 
discord might have a role to play. We will not delve into a general analysis of this; 
rather we will point out the role of quantum discord through a simple example. We 
also thank W. H. Zurek for pointing out the problem. 

The phenomenon we have in mind is the distribution of entanglement. The 
significance of this task cannot be overemphasized in quantum information science. 
To be able to entangle two distant particles is at the heart of quantum communication 
and computation. That this can be achieved via a third mediating particle which 
remains separable at all times from the other two was shown in [ ] . Let us, in very 
brief, demonstrate this claim. 

We start with a tripartite state pabc, where a and b denote the systems of Alice 
and Bob, while c is the mediator particle or channel, which starts out in Alice's 
possession. In our particular example, taken from [in]. 



A:=0 i=0 



where 



|0>+e"""^|l) 



This state is separable between all three parties as it is a convex combination of 
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Pabc 



(2.47) 



product states. In the standard product basis, 

/ioooooio\ 
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1 00000000 
6 00001000 
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Firstly, Alice performs a controlled NOT (CNOT) gate operation between a and 
c (with a as the control), producing the state 



(^abc 



(2.48) 



1 1 

-(|000) + |111))((000| + (111|) + - J2 P'.3k\tjk){tjk\, 

i,j,k=0 

with all /3's being except for /9ooi = /^oio = Aoi = Aio = 1- The particle c is then 
sent over to Bob. Since the CNOT was executed between a and c, the state must 
be separable with respect to the b-ac split, cxab and a^c are also separable. Also, the 
state is separable under the c-ab split as the state is invariant under the exchange of 
b and c. Likewise, (Tac is separable. The a-bc split is more interesting though. For 
this we write in the standard product basis, 

/lOOOOOOl\ 
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1 00000000 
6 00000000 

10 
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10 1 



O'abc 



(2.49) 
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Then the spectrum of the partial transpose of (Xabc with respect to a is 



11111111 



6 5 /^' /^' /^' /^'/^' /^' /"' 




which has a negative eigenvalue. Hence, the a-bc split actually contains entanglement. 

Once Bob receives the particle, he does a CNOT between b and c, with the former 
as the control. The resulting state is 



1 



'^abc 



$H 



f$+|®|o)(o| + -u®|i)(i|. 



(2.50) 



3 ' ' ^ ' ' ' ^ ' 3 
where |$+) = (|00) + \\\))/\/2 is the maximally entangled state. In the standard 
product basis, 

/ioooooio\ 
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(2.51) 
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The particle c is still separable from the other two particles, which we have managed 
to get entangled. However, the a-bc split and the b-ac split have the same amount 
of entanglement as in the a-bc split of (Xabc- One uses this entanglement to get a Bell 
state by measuring c in the standard basis. On average, we can thus use this 1 ebit 
of entanglement 1/3 of the time for any desired purpose. 

Having demonstrated our claim, we now invest in a little introspection. If entan- 
glement is to behave as a resource, one would expect some nature of conservation 
law to hold. At least, we should be able to argue that the entanglement we gen- 
erated was a result of some form of expenditure. We performed two CNOT gates, 
which generates some entanglement in the intermediate states. However, we are at 
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a loss to take this any further due to a lack of our ability quantifying entanglement, 
particularly in entropy units. We will now show that quantum discord will provide 
a natural accounting of the resources and a more soothing resolution. 

Let us consider measurements on c, as that is the most contentious party in our 
protocol. The projectors are 

Hi = lab ® ^^^ , n2 = lab ® -^^ , (2.52) 

with tti = sin 6 cos 0, 02 = sin ^ sin 0, 0:3 = cos 6 as usual. 

We will start with the state pabc- A simple calculation shows 

3 + cos 9 3 — cos 9 

P^ = a ' ^2 = 7 , 

o 

and that the spectrum of the reduced operators are 

, , 1 , f 11 + cos 9 1 + cos 9 1 — cos 9 ^ , ^ 

\{pl,) = 0,0,0,0,-,—— ^,;^77— ^,^^777- ^ , 2.53 

[ 2 2(3 + cos t^j 2(3 + cos t^j 2(3 + cos t^j J 

2, r^^^^l 1 + cos 9 1 — cos 9 1 — cos I 



\{pi,) = 0,0,0,0,-,— -,— -,— . (2.54) 

[ 2 2(3 — cos t^j 2(3 — cos t^j 2(3 — cos t^) J 

We immediately see that the minimization in this case is a single variable affair (0 
having dropped out), and this quickly leads to (for ^ = 0) 

V{pabc) = 0. (2.55) 

Next, let us calculate the discord for the state (Jabc- We find 

1 

Pi = P2 = 2' 

and that the spectrum of the reduced operators is 

,,1, ,,0, f 1 o^ 1 <,9 2 + sin^ 2-sin^l 

\{al,) = \{al,) = |0, 0, 0, 0, - cos^ -, - sin^ -, ^ , ^ j . (2.56) 

Again, the minimization involves involves one variable, whose minimum is attained 
for ^ = 0, as 

V{aabc) = \- (2.57) 
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Finally, a similar calculation for Tabc provides 

-DiTabc) = 0. (2.58) 

The three results, in Eqs 2.55, 2.57, 2.58 lead us to very striking conclusions, viz., 

1. The amount of discord generated is exactly equal to the average amount of 
entanglement we can extract in the next stage of the protocol. 

2. The fact there is no discord at the end of the protocol is circumstantial evidence 
that the discord is somehow 'converted' into entanglement. 

This might be construed to imply that discord is a more fundamental resource for 
quantum information than entanglement. In our view, however, further investiga- 
tion is called for before such conclusions are drawn. Indeed, a genuine measure of 
tripartite entanglement may help clear up the resource nature of entanglement in 
the protocol. Another line of investigation would be to demonstrate the qualitative 
and/or quantitative role of discord in the general scheme of entanglement distribu- 
tion [ ], of which we presented just an example, or in other mixed-state protocols 
involving entanglement. A third line of research would be to arrive at an independent 
operational interpretational for quantum discord. An interpretation of this nature 
will demonstrate the power, utility and resourcefulness of quantum discord. 

With these statements, we will move on to our study of mixed-state quantum 
computation in Chapters 3 and 4. 
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Chapter 3 



Entanglement in the DQCl model 



Why then, can one desire too much of a good thing ? 



-William Shakespeare in As You Like It 



The purpose of this chapter is to investigate the existence of and amount of 
entanglement in the DQCl circuit that is used to estimate the normalized trace. The 
DQCl model consists of a special qubit (qubit 0) in the initial state |0)(0| = |(/i + Z), 
where Z is a Pauli operator, along with n other qubits in the completely mixed state, 
J„/2", which we call the unpolarized qubits. The circuit consists of a Hadamard gate 
on the special qubit followed by a controlled unitary on the remaining qubits [57]: 



|0)(0| ^H 

In/2'' 



(^= ^" 



Figure 3.1: The DQCl circuit 
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After these operations, the state of the n + 1 qubits becomes 

Pn+l = ^(|0)(0|®J„+|l)(l|®J„+|0)(l|®f/t+|l)(0|®f/„) = ^ I ^" ^" 

(3.1) 

where A^ = 2". The information about the normahzed trace of f/„ is encoded in 
the expectation values of the Pauh operators X and Y of the special qubit, i.e., 

(X) = Re[tr(f/„)]/2" and {¥) = -Im[tr(f/„)]/2". 

To read out the desired information, say, about the real part of the normalized 
trace, one runs the circuit repeatedly, each time measuring X on the special qubit 
at the output. The measurement results are drawn from a distribution whose mean 
is the real part of the normalized trace and whose variance is bounded above by 
1. After L runs, one can estimate the real part of the normalized trace with an 
accuracy e ~ Xj^jL. Thus, to achieve accuracy e requires that the circuit be run 
L ~ 1/e^ times. More precisely, what we mean by estimating with fixed accuracy 
is the following: let Pe be the probability that the estimate is farther from the true 
value than e ; then the required number of runs is L ~ ln(l/Pe)/e^- That the number 
of runs required to achieve a fixed accuracy does not scale with number of qubits 
and scales logarithmically with the error probability is what is meant by saying that 
the DQCl circuit provides an efficient method for estimating the normalized trace. 

Throughout much of our analysis, we use a generalization of the DQCl circuit, 
in which the initial pure state of the special qubit is replaced by the mixed state 
\{l\ + olZ\ which has polarization a, giving an overall initial state 






U. 



n 



Figure 3.2: The DQCl circuit with arbitrary polarization 
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(3.2) 



We generally assume that a > 0, except where we explicitly note otherwise. After 
the circuit is run, the system state becomes 
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(3.3) 



The effect of subunity polarization is to reduce the expectation values of {X) and (Y) 
by a factor of a, thereby making it more difficult to estimate the normalized trace. 
Specifically, the number of runs required to estimate the normalized trace becomes 
L ~ \n{l/Pe)/a^e^. Reduced polarization introduces an additional overhead, but 
as long as the special qubit has nonzero polarization, the model still provides an 
efficient estimation of the normalized trace. What we are dealing with is really the 
"power of even the tiniest fraction of a qubit." 

For n+l qubits, all states contained in a ball of radius r„+i centered at the 
completely mixed state are separable [9, 39] (distance is measured by the Hilbert- 
Schmidt norm). Unitary evolution leaves the distance from the completely mixed 
state fixed, so at all times during the circuit in Fig. (3.2), the system state is a fixed 
distance A/ti^Pi— Xm/2iV)^ = a2~*^"+^)/^ from the completely mixed state. This 
suggests that with a small enough, there might be an exponential speedup with 
demonstrably separable states. This suggestion doesn't pan out, however, because 
the radius of the separable ball decreases exponentially faster than 2^*^"^^-*'^. The 
best known lower bound on r„ is 2 x 6~"/^ [ ] ; for the system state to be contained 
in a ball given by this lower bound, we need a < 2 x 3^*^"'^^^^^. The exponential 
decrease of a means that an exponentially increasing number of runs is required to 
estimate the normalized trace with fixed accuracy. More to the point, the possibility 
that the actual radius of the separable ball might decrease slowly enough to avoid 
an exponential number of runs is ruled out by the existence of a family of n-qubit 
entangled states found by Diir et al. [ ], which establishes an upper bound on r„ 
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that goes as 2 x 2~"' for large n, implying that a < 2 x 2^'-""^^)/^ if the system state is 
to be in the ball given by the upper bound. These considerations do not demonstrate 
the impossibility of an exponential speedup using separable states, but they do rule 
out the possibility of finding such a speedup within the maximal separable ball about 
the completely mixed state. 

We are thus motivated to look for entanglement in states of the form in Eq. (3.3), 
for at least some unitary operators f/„. Initial efforts in this direction are not en- 
couraging. It is clear from the start that the marginal state of the n unpolarized 
qubits remains completely mixed, so these qubits are not entangled among them- 
selves. Moreover, in the state in Eq. (3.3), as was shown in Ref. [(i(i], the special 
qubit is unentangled with the n unpolarized qubits, no matter what Un is used. To 
see this, one plugs the eigendecomposition of the unitary, Un = J2j ^^'^^\^j){^j\y i^^o 
the expression for p„+i(a). This gives a separable decomposition 

Pn+i(a) = ^ '^{\aj){aj\ + \bj){bj\) ® \ej){ej\ , (3.4) 

j 

where \aj) = cos6'|0) + e'^^^ sin6'|l) and \bj) = sin6'|0) + e'^^^ cos 6^1 1), with sin26' = a. 

No entanglement of the special qubit with the rest and no entanglement among 
the rest — where then are we to find any entanglement? We look for entanglement 
relative to other divisions of the qubits into two parts. In such bipartite divisions the 
special qubit is grouped with a subset of the unpolarized qubits. To detect entan- 
glement between the two parts, we use the Peres-Horodecki partial transpose crite- 
rion [64, 45] , and we quantify whatever entanglement we find using a closely related 
entanglement monotone which we call the multiplicative negativity [ ]. The Peres- 
Horodecki criterion and the multiplicative negativity do not reveal all entanglement — 
they can miss what is called bound entanglement — but we are nonetheless able to 
demonstrate the existence of entanglement in states of the form in Eq. (3.1) and 
Eq. (3.3). For convenience, we generally refer to the multiplicative negativity sim- 
ply as the negativity. The reader should note, as we discuss in Sec. 2.1, that the 
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term "negativity" was originally applied to an entanglement measure that is closely 
related to, but different from the multiplicative negativity. 



The amount of entanglement depends, of course, on the unitary operator f/„ and 
on the bipartite division. We present three results in this regard. First, in Sec. 3.2, 
we construct a family of unitaries Un such that for a > 1/2, pn+i{<^) is entangled for 
all bipartite divisions that put the first and last unpolarized qubits in different parts, 
and we show that for all such divisions, the negativity is (2a + 3)/4 for a > 1/2 (5/4 
for a = 1), independent of n. Second, in Sec. 3.3, we present numerical evidence that 
the state p„+i of Eq. (3.1) is entangled for typical unitaries, i.e., those created by 
random quantum circuits. For n + 1 = 5, . . . , 10, we find average negativities between 
1.155 and just above 1.16 for the splitting that puts \n/2\ of the unpolarized qubits 
with qubit 0. Third, in Sec. 3.4, we show that for all unitaries and all bipartite 
divisions of the n + 1 qubits, the negativity of p„+i(a) is bounded above by the 
constant \^1 + a^ {\f2 ^ 1.414 for a = 1), independent of n. Thus, when n is large, 
the negativity achievable by the DQCl circuit (3.1) becomes a vanishingly small 
fraction of the maximum negativity, ~ 2"/^, for roughly equal bipartite divisions. 



The layout of this chapter is as follows. In Sec. 3.1 we examine the classical 
problem of estimating the normalized trace of a unitary. In Sec. 2.1 we reviewed the 
pertinent properties of the negativity. We apply these to obtain our three key results 
in Sees. 3.2-3.4. We conclude in Sec. 3.5 and prove a brief Lemma in Appendix B. 
Throughout we use A to stand for the partial transpose of an operator A relative to 
a particular bipartite tensor-product structure, and we rely on context to make clear 
which bipartite division we are using at any particular point in this chapter. 



Chapter 3. Entanglement in the DQCl model 51 

3.1 Classical evaluation of the trace 

Every fair judgement on the power of quantum computation to solve a mathematical 
problem must be made relative to the best known classical algorithm to solve the 
same. In that spirit, this section briefly outlines a classical method for estimating the 
trace of a unitary operator that can be implemented efficiently in terms of quantum 
gates, and we indicate why this appears to be a problem that is exponentially hard 
in the number of qubits. 

The trace of a unitary matrix f/„ = f/ is the sum over the diagonal matrix 
elements of U: 



(f/) = 5^(a|f/|a) . (3.5) 



tr 



Here a is a bit string that specifies a computational-basis state of the n qubits. By 
factoring U into a product of elementary gates from a universal set and inserting a 
resolution of the identity between all the gates, we can write tr(f/) as a sum over the 
amplitudes of Feynman paths. A difficulty with this approach is that the sum must 
be restricted to paths that begin and end in the same state. We can circumvent this 
difficulty by preceding and succeeding U with a Hadamard gate on all the qubits. 
This does not change the trace, but does allow us to write it as 

tr(f/) = J2 (a| ^®" |b) (b| U |c) (c| i/®" \^) = \Y1 (-1)"^^^'^ (b| U |c) . (3.6) 

a,b,c a,b,c 

Now if we insert a resolution of the identity between the elementary gates, we get 
tr(f/) written as an unrestricted sum over Feynman-path amplitudes, with an extra 
phase that depends on the initial and final states. 

Following Dawson et al. [' *^*], we consider two universal gate sets: (i) the Hadamard 
gate H, the 7r/4 gate T, and the controlled-NOT gate and (ii) H and the Toffoli gate. 
With either of these gate sets, most of the Feynman paths have zero amplitude. Daw- 
son et al. [19] introduced a convenient method, which we describe briefly now, for 



Chapter 3. Entanglement in the DQCl model 52 

including only those paths with nonzero amplitude. One associates with each wire 
in the quantum circuit a classical bit value corresponding to a computational basis 
state. The effect of an elementary gate is to change, deterministically or stochasti- 
cally, the bit values at its input and to introduce a multiplicative amplitude. The 
two-qubit controUed-NOT gate changes the input control bit x and target bit y de- 
terministically to output values x and y(Bx, while introducing only unit amplitudes. 
Similarly, the three-qubit Toffoli gates changes the input control bits x and y and 
target bit z deterministically to x, y, and z © xy, while introducing only unit am- 
plitudes. The T gate leaves the input bit value x unchanged and introduces a phase 
gia;7r/4_ rjj^^ Hadamard gate changes the input bit value x stochastically to an output 
value y and introduces an amplitude {—l)^^/^/2. 

The classical bit values trace out the allowed Feynman paths, and the product of 
the amplitudes introduced at the gates gives the overall amplitude of the path. In 
our application of evaluating the trace (3.6), a path is specified by n input bit values 
(which are identical to the output bit values), n random bit values introduced by 
the initial Hadamard gates, and h random bit values introduced at the h Hadamard 
gates required for the implementation of U. This gives a total oi 2n + h bits to 
specify a path and thus 2^""'"^ allowed paths. We let x denote collectively the 2n + h 
path bits. 

If we apply the gate rules to a Hadamard- Toffoli circuit, the only gate amplitudes 
we have to worry about are the ±l/-\/2 amplitudes introduced at the Hadamard 
gates. There being no complex amplitudes, the trace cannot be complex. Indeed, 
for this reason, achieving universality with the if -Toffoli gate set requires the use 
of a simple encoding, and we assume for the purposes of our discussion that this 
encoding has already been taken into account. With all this in mind, we can write 
the trace (3.6) as a sum over the allowed paths, 

Mf/) = ^E(-1)'^^^- (3-7) 

X 
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Here ip^x.) is a polynomial over Z2, specifically, the mod-2 sum of the products of 
input and output bit values at each of the Hadamard gates. The downside is that 
a string of TofFoli gates followed by a Hadamard can lead to a polynomial that is 
high order in the bit values. As pointed out by Dawson et al. [ ], we can deal with 
this problem partially by putting a pair of Hadamards on the target qubit after each 
Toffoli gate, thus replacing the quadratic term in the output target bit with two new 
random variables and preventing the quadratic term from iterating to higher order 
terms in subsequent Toffoli gates. In doing so, we are left with a cubic term in V^(x) 
from the amplitude of the first Hadamard. The upshot is that we can always make 
ipi'x) a cubic polynomial. 



tr(f/) ^ 



2n+h/2 



(3.8) 



Notice now that we can rewrite the trace as 

number of x such \ / number of x such 
that V(x) =0 J \ that ^{x.) = 1 

thus reducing the problem of evaluating the trace exactly to counting the number of 
zeroes of the cubic polynomial ipi'x). This is a standard problem from computational 
algebraic geometry, and it is known that counting the number of zeroes of a general 
cubic polynomial over any finite field is #P complete [25]. It is possible that the 
polynomials that arise from quantum circuits have some special structure that can 
be exploited to give an efficient algorithm for counting the number of zeroes, but in 
the absence of such structure, there is no efficient classical algorithm for computing 
the trace exactly unless the classical complexity hierarchy collapses and all problems 
in #P are efficiently solvable on a classical computer. 

Of course, it is not our goal to compute the trace exactly, since the quantum 
circuit only provides an efficient method for estimating the normalized trace to fixed 
accuracy. This suggests that we should estimate the normalized trace by sampling 
the amplitudes of the allowed Feynman paths. The normalized trace. 



ti{U) _ 

On 9271 



l_^2V^(_l),(x)^ (3.9) 
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which hes between —1 and +1, can be regarded as the average of 2^"'"^^ quantities 
whose magnitude, 2^^/^, is exponentially large in the number of Hadamard gates. To 
estimate the average with fixed accuracy requires a number of samples that goes as 
2^, implying that this is not an efficient method for estimating the normalized trace. 
The reason the method is not efficient is pure quantum mechanics, i.e., that the trace 
is a sum of amplitudes, not probabilities. 

If we apply the gate rules to a Hadamard- T-controUed-NOT circuit, the bit value 
on each wire in the circuit is a mod-2 sum of appropriate bit values in x, but now 
we have to worry about the amplitudes introduced by the Hadamard and T gates. 
The trace (3.6) can be written as 

tr(f/) = ^^^^^^^'^"^''^(-1)'^''^ ■ (3.10) 

X 

Here 0(x) is a polynomial over Z2, obtained as the mod-2 sum of the products of 
input and output bit values at each of the Hadamard gates. Since the output value 
is a fresh binary variable and the input value is a mod-2 sum of bit values in x, (/>(x) 
is a purely quadratic polynomial over Z2. The function x(x) is a mod-8 sum of the 
input bit values to all of the T gates. Since these input bit values are mod-2 sums 
of bit values in x, x(x) is linear in bit values, but with an unfortunate mixture of 
mod-2 and mod-8 addition. We can get rid of this mixture by preceding each T gate 
with a pair of Hadamards, thus making the input to the every T gate a fresh binary 
variable. With this choice, x(x) becomes a mod-8 sum of appropriate bit values 
from X. 



We can rewrite the sum (3.10) in the following way 



2n+h/2 

j=0 



(,^JJ 



number of x such that \ I number of x such that 
X(x) = j and 0(x) = Oj \x{^) = 3 and 0(x) = 1^ 

(3^11) 



Thus the problem now reduces to finding simultaneous (binary) solutions to the 
purely quadratic Z2 polynomial 0(x) and the purely linear Zg polynomial x(x). One 
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has to be careful here to note that we are only interested in binary solutions, so we 
are not solving x(x) = J over all values in Zg. The number of solutions of a purely 
quadratic polynomial over Z2 can be obtained trivially [ ], but the constraint over 
Zg means that one must count the number of solutions over a mixture of a field 
and a ring. The complexity class for this problem is not known, but given the 
equivalence to counting the number of solutions of a cubic polynomial over Z2, it 
seems unlikely that there is an efficient classical algorithm. Moreover, an attempt to 
estimate the normalized trace by sampling allowed paths obviously suffers from the 
problem already identified above. 



3.2 Entanglement in the DQCl circuit 

In this section we present our main results. As outlined earlier, we will use the 
negativity as measure of entanglement. The negativity is the sum of the singular 
values of p. For states of the form we are interested in, given by Eq. (3.3), the 
negativity is determined by the singular values of the partial transpose of the unitary 
operator f/„. To see this, consider any bipartite division of the qubits. Performing 
the partial transpose on the part that does not include the special qubit, we have 



Pn+l(«) = ^ " " h (3-12) 




where f7„ is the partial transpose of f/„. relative to the chosen bipartite division. 
Notice that if we make our division between the special qubit and all the rest, then 
Un = f/J is a unitary operator, and p„+i(a) is the quantum state corresponding to 
using f/J in the circuit (3.2); this shows that for this division, the negativity is 1, 
consistent with our earlier conclusion that the special qubit is not entangled with the 
other qubits. For a general division, we know there are unitaries V and W such that 
Un = VSW, where S is the diagonal matrix of singular values Sj{Un)- This allows 



Chapter 3. Entanglement in the DQCl model 56 

us to write 

Pn+i(a) = \^\ J ' (313) 

showing that p„+i(a) is a unitary transformation away from the matrix in the middle 
and thus has the same eigenvalues. The block structure of the middle matrix makes 
it easy to find these eigenvalues, which are given by [1 ± aSj(f7„)]/2A^. This allows 
us to put the negativity in the form 

1 ^ . . 1 ^ 

M{pn+iia)) = —^\l + aSj{Un)\ + \l-aSj{Un)\ = — ^max(|a| sj(f/„), l) , 

(3.14) 






which is valid for both positive and negative values of a. An immediate consequence 
of Eq. (3.14) is that Ai(^pn+i{a)) = A^(p„+i(— a)), as one would expect. Since 
Pn+i(«) is a mixture of p„_|_i(+l) = Pn+i and p„+i(— 1), convexity tells us immediately 
that A^(p„+i(a)) < Ai[pn+i), i.e., that a mixed input for the special qubit cannot 
increase the negativity over that for a pure input. More generally, we have that the 
negativity cannot decrease at any point as a increases from to 1. 

As the first result of this section, we construct a family of unitaries Un that 
produce global entanglement in the DQCl circuit (3.1). For a = 1, the negativity 
produced by this family is equal to 5/4, independent of n, for all bipartite divisions 
that put the first and last unpolarized qubits in different parts. We conjecture that 
this is the maximum negativity that can be achieved in a circuit of the form (3.1). 

Before the measurement, the output state of the circuit (3.2) is given by Eq. (3.3). 
To construct the unitaries f/„, we first introduce a two-qubit unitary matrix 



U2^\ , (3.15) 




where Ai, Bi, Ci, and Di are single-qubit (2 x 2) matrices that must satisfy A\Ai + 
dIDi = bIBi + ClCi = h and A\Ci + D\Bi = to ensure that U2 is unitary. The 
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n-qubit unitary f/„ is then defined by 

'J n = I 

+ |0)(1| ® X„_2 ® Ci + |1)(0| ® Xn-2 ® 1^1 . (3.16) 

Here we use Xi, Yi, and Zi to denote single-qubit Pauli operators. A subscript k 
on the identity operator or a Pauh operator denotes a tensor product in which that 
operator acts on each of k qubits. If we adopt the convention that Xq = Jq = 1, then 
Un reduces to f/2 when n = 2. It is easy to design a quantum circuit that realizes [/„■ 
The structure of the circuit is illustrated by the case of f/4: In general, the two-qubit 

1st qubit 



4th qubit 

2nd qubit -^ 



Uo 



-^ 



3rd qubit © ® 

Figure 3.3: Unitary t/4 that generates entanglement in DQCl circuit for all n 

unitary U2, acting on the first and last qubits, is bracketed by controlled-NOT gates 
from the first qubit, acting as control, to each of the other qubits, except the last, 
as targets. 

Because /i and Xi are invariant under transposition, it is clear from the form 
of Un that in the state in Eq. (3.3), all qubits, except 0, 1, and n, are invariant 
under transposition. We can use this fact to find the negativity for all bipartite 
divisions. First consider any bipartite division that puts qubits 1 and n in the same 
part. There are two possibilities. If the special qubit is in the same part as qubits 1 
and n, then partial transposition on the other part leaves p„+i(a;) unchanged, so the 
negativity is 1. If the special qubit is not in the same part as 1 and n, then partial 
transposition on the part that includes 1 and n is the same as partial transposition 
of all the unpolarized qubits, a case we already know to have negativity equal to 
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1. We conclude that any bipartite division that puts 1 and n in the same part has 
negativity equal to 1. 

Turn now to bipartite divisions that put qubits 1 and n in different parts. There 
are two cases to consider: (i) the special qubit is in the same part as qubit 1, and 
(ii) the special qubit is in the same part as qubit n. In case (i), partial transposition 
of the part that contains qubit n gives 

Pn+l(«) = 7^r\^ With Un= \ \ ■ (3.17) 

2^ \aUn In J \X„„2 ® Dj In-2 ® Bf 

The structure of f/„ comes about due to the transposition over the last qubit in 
Eq. (3.16). In case (ii), partial transposition of the part that contains qubit 1 gives 
(here the transposition is over the first qubit in Eq. (3.16)) 

\Xn^2^Ci In-2®Bi 

The basic structure of p„+i(q;) is the same in both cases. Without changing the 
spectrum, we can reorder the rows and columns to block diagonalize Pn+iio) so that 
there are N /A blocks, each of which has the form 

^P3(«) , (3.19) 

where pz{a) is the appropriate partial transpose of the three-qubit output state. 
Thus the spectrum of p„+i(a) is the same as the spectrum of P3(a), except that 
each eigenvalue is reduced by a factor of A/N . In calculating the negativity, since 
each eigenvalue is (A^/4)-fold degenerate, the reduction factor of 4/A^ is cancelled by a 
degeneracy factor of A^/4, leaving us with the fundamental result of our construction. 



M[pn+i{a)) = M[p^{a)) . (3.20) 
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This applies to both cases of bipartite sphttings that we are considering, showing 
that all divisions have the same negativity as the corresponding n = 2 construction. 

We now specialize to a particular choice of f/2 given by 



A, 




Bi 




Ci 




and Di 




For this choice, the two cases of bipartite division lead to the same partial transpose. 
The spectrum of 

f h aAi aDi\ 



P3[a) 



Ji aCi aBi 
aAi aDi Ii 
\aCi aBi h J 



(3.22) 



IS 



Spec(p3(a)) = -(1 + 2a, 1, 1, 1, 1, 1, 1 - 2a) , 
giving a negativity equal to 1 for a < 1/2 and a negativity 
M{pn+i{a)) =M{p3{a)) = -{2a + 3) for a > 1/2 . 



(3.23) 



(3.24) 



This result shows definitively that the circuit can produce entanglement, at least 
for a > 1/2. We stress that the negativity achieved by this family of unitaries is 
independent of n > 2. 

For a = 1, the negativity achieved by this family reduces to 5/4. For large n, 
this amount of negativity is a vanishingly small fraction of the maximum possible 
negativity, ~ 2"/^, for roughly equal divisions of the qubits. This raises the question 
whether it is possible for other unitaries to achieve larger negativities. A first idea 
might be to find two-qubit unitaries f/2 that yield a higher negativity A4{p3) = 
A4{pn+i) when plugged into the construction of this section, but the bounds we 
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find in Sec. 3.4 dispose of tliis notion, since tfiey sliow tliat 5/4 is tlie maximum 
negativity that can be acliieved for n = 2. Anotlier approach would be to generahze 
the construction of this section in a way that is obvious from the circuit in Fig. (3.3), 
i.e., by starting with a fc-qubit unitary in place of the two-qubit unitary of Fig. (3.3). 
Numerical investigation of the case k = 3 has not turned up negativities larger than 
5/4. We conjecture that 5/4 is the maximum negativity that can be achieved by 
states of the form (3.1). Though we have not been able to prove this conjecture, we 
show in the next section that typical unitaries for n + 1 < 10 achieve negativities 
less than 5/4 and in the following section that the negativity is rigorously bounded 
by v^. 

We stress that we are not suggesting that the construction of this section, with 
U2 given by Eq. (3.21), achieves the maximum negativity for all values of a, for 
that would mean that we believed that the negativity cannot exceed 1 for a < 1/2, 
which we do not. Although we have not found entanglement for a < 1/2, we suspect 
there are states with negativity greater than 1 as long as a is large enough that 
Pn+i{a) lies outside the separable ball around the maximally mixed state, i.e., a > 



2("+i)/V„,+i. The bound of Sec. 3.4 only says that M{pn+iia)) < Vl + a^ thus 
allowing negativities greater than 1 for all values of a except a = 0. Moreover, since 
the negativity does not detect bound entanglement, there could be entangled states 
that have a negativity equal to 1. 

3.3 The average negativity of a random unitary 

Having constructed a family of unitaries that yields a DQCl state with negativity 
5/4, a natural question to ask is, 'What is the negativity of a typical state produced 
by the circuit in Fig. (3.1)?" To address this question, we choose the unitary oper- 
ator in the circuit in Fig. (3.1) at random and calculate the negativity. Of course, 
one must first define what it means for a unitary to be "typical" or "chosen at 
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random". The natural measure for defining this is the Haar measure, which is the 
unique left-invariant measure for the group U(A^) [13]. The resulting ensemble of 
unitaries is known as the Circular Unitary Ensemble, or CUE, and it is parameter- 
ized by the Hurwitz decomposition [68]. Although this is an exact parameterization, 
implementing it requires computational resources that grow exponentially in the size 
of the unitary [ ] . To circumvent this, a pseudo-random distribution that requires 
resources growing polynomially in the size of the unitary was formulated and investi- 
gated in Ref. [')(']. This is the distribution from which we draw our random unitaries, 
and we summarize the procedure for completeness. 

We first define a random SU(2) unitary as 

(e*"^ cos 6 e'^ sin 6 \ 
(3.25) 
-e-*^sin^ e-'^ cose J 

where 6 is chosen uniformly between and 7r/2, and and x ^^^ chosen uniformly 
between and 27r. A random unitary applied to each of the n qubits is then 

n 

R = <^R{e„<Pi,Xi), (3.26) 

where a separate random number is generated for each variable at each value of i. 
Now define a mixing operator M in terms of nearest-neighbor Z ^ Z couplings as 

n-l 
4 



M = exp I z^ ^ Z(^) ® Z(J+i) ) . (3.27) 



The pseudo-random unitary is then given by 

RjMRj^i-.-MRaMRi , (3.28) 

where j is a positive integer that depends on n, and each R^ is chosen randomly 
as described above. For a given n, the larger j is, the more accurately the pseudo- 
random unitary distribution resembles the actual CUE. From the results in Ref. [30] , 
j = 40 gives excellent agreement with the CUE for unitary operators on at least up 
to 10 qubits, so this is what we use in our calculations. 
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Average Negativity v. Number olQubits 
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Figure 3.4: Average negativity of the state p„_|_i of Eq. (3.1) [a = 1) for a randomly 
chosen unitary f/„ for two different bipartite sphttings, (n, 1) and {[n/2\ + 1, \n/2]). 
The (n, 1) sphtting appears to reach an upper bound quickly, whereas the other 
sphtting is still rising slowly at 10 qubits. 



Due to boundary effects, not all bipartite splittings that put k unpolarized qubits 
in one part are equivalent. Nevertheless, we consider only bipartite divisions that 
split the qubits along horizontal lines placed at various points in the circuit of 
Eq. (3.1). We refer to the division that groups the last k qubits together as the 
{n + 1 — k,k) splitting. For a = 1, we calculate the average negativity and standard 
deviation of a pseudo-random state pn+i for two different bipartite splittings, (n, 1) 
and {[n/2\ + 1, \n/2]). These results are plotted in Figs. 3.4 and 3.5. 

For n + 1 = 5, . . . , 10, the average negativity for the roughly equal splitting lies 
between 1.135 and just above 1.16. The standard deviation appears to converge 
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Standard Deviation v. Number of Qubits 
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Figure 3.5: Semi-log plot of the standard deviation in the negativity of the randomly 
chosen state pn+i- The fit curves show that the standard deviation is decaying 
exponentially, so that for large numbers of qubits, almost all unitaries give the same 
negativity. 



exponentially to zero, as in Ref. [72], a behavior that is typical of asymptotically 
equivalent matrices. In addition, for 9 + 1 qubits, we calculate the average negativity 
and standard deviation for all nontrivial {k ^ n) bipartite splittings {n + 1 — k,k), 
and the results are shown in Fig. 3.6. 



3.4 Bounds on the negativity 



In this section, we return to allowing the special qubit in the circuit (3.2) to have 
initial polarization a. Since the value of n is either clear from context or fixed, we 
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Figure 3.6: Average negativity of the state pio of Eq. (3.1) (a = 1) for a randomly 
chosen unitary Ug as a function of bipartite sphtting number k for bipartite sphttings 
(10 — k,k). The error bars give the standard deviations. The function attains a 
maximum when the bipartite spht is made between half of the qubits on which the 
unitary acts. 



reduce the notational clutter by denoting the state p„+i(a) of Eq. (3.3) as pa- 

Given a particular bipartite division, the partial transpose of pa with respect to 
the part that does not include the special qubit is 

° "■- (3.29) 




Using the binomial theorem, we can expand tr(p^) in terms of tr(C ): 

tr(p^)=(^yE(l)«'t^(^')- (3-30) 

^ ^ fc=0 ^ ^ 
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When k is odd, C'' is block off-diagonal, so its trace vanishes. When k is even, we 
have 

tr(C"=) = 2tr((f/„[/t)'=/2). (3.31) 

When k = 2, this simplifies to tT{C^) = 2tr(f/„[/t) = 2tr(t/„f/t) = 2N. The crucial 
step here follows immediately from the property tT{AB) = tT{AB), which we prove 
as a Lemma in Appendix B. Note that in general tr(y4iy42 . . . Ai) ^ tT{AiA2 . . . Ai) if 
/ > 2, so we cannot give a similar general calculation of tr(p^) for even s > 4, since 
it involves terms of this form. 

Using Eq. (3.30), we can now obtain three independent constraint equations on 
the eigenvalues \j = \j{pa) of the partial transpose pa'- 

2N ^ 

Y^ A • = tr(p^) = ^^^^^[(1 + ay + (1 - ay] , . = 1, 2, 3. (3.32) 

Since the negativity is given by 

M{pa) = J2\^^\ ' (3.33) 

i 

we can find an upper bound on the negativity by maximizing ^ . |Aj| subject to the 
constraints (3.32). If we consider only the s = 1, 2 constraints, we obtain a nontrivial 
upper bound on the negativity with little effort. We find that adding the constraint 
s = 3 adds nothing asymptotically for large N, but for small A^ yields a tighter 
bound than we get from the s = 1, 2 constraints, although this comes at the cost 
of considerably more effort. We emphasize that these bounds apply to all bipartite 
divisions and to all unitaries [/„. Notice that we have no reason to expect these 
bounds to be saturated, since the traces of higher powers of pa impose additional 
constraints that we are ignoring. The one exception is the case of three qubits, where 
the s = 1, 2, 3 constraints are a complete set, and indeed, in this case, the s = 1, 2, 3 
bound is 5/4, which is saturated by the unitary found in Sec. 3.2. 

The remainder of this section is devoted to calculating the s = 1,2 and s = 1, 2, 3 
upper bounds. A graphical summary of our results for a = 1 is presented in Fig. 3.7. 
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Figure 3.7: Plot of the bounds on the negativity of states of the form (3.1), i.e., for 
a pure-state input in the zeroth register (a = 1). The uppermost plot is the simple 
analytic bound A1i,2 = "\/2, obtained using the s = 1, 2 constraint equations; the 
next largest plot is the numerically constructed s = 1, 2, 3 bound. One can see that 
the s = 1, 2, 3 bound asymptotes to the s = 1, 2 bound. As noted in the text, these 
bounds are independent of the unitary [/„ and the bipartite division. The flat line 
shows the negativity 5/4 for the state constructed in Sec. 3.2, currently the state of 
the form (3.1) with the largest demonstrated negativity; notice that for n + 1 = 3, 
this state attains the s = 1, 2, 3 bound. The lowest two sets of data points display 
the expected negativities for a randomly chosen unitary using the bipartite splittings 
(n, 1) and ([n/2j + 1, \n/2\), which were also plotted in Figs. 3.4 and 3.5. 
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3.4.1 The s = 1,2 Bound 

We can use Lagrange multipliers to reduce the problem to maximizing a function of 
one variable, but first we must deal with the absolute value in Eq. (3.33). To do so, we 
assume that t of the eigenvalues are negative and the 2N — t others are nonnegative, 
where t becomes a parameter that must now be included in the maximization. We 
want to maximize 

t 2N 

-Mi,2 = -J]A,+ Y, ^^ (3-34) 

i=i j=t+i 

subject to the constraints 

2N ^^ 1 I 2 

J2X, = 1 and E^^ = ^- (3.35) 

fc=i fc=i 

The notation we adopt here for the indices is that i labels negative eigenvalues and 
j labels nonnegative eigenvalues, while k can label either. This serves to remind us 
of the sign of an eigenvalue just by looking at its index. 

Introducing Lagrange multipliers fi and u, the function we want to maximize is 

t 2N / 2N \ / 2N 1 I 2 \ 

/(A.,t) = -5^A.+ E^^+Z^ E^'^-1 +'^ E^^-^ • (3-36) 

j=l j=t+l \fc=l / \A:=1 / 

Differentiating with respect to Aj and then Xj, we find 

- 1 + /i + 2uXi = , (3.37) 

+1 + /i + 2uXj = . (3.38) 

We immediately see that in the maximal solution, all the negative eigenvalues are 
equal, and all the nonnegative eigenvalues are equal. We can now reformulate the 
problem in the following way. If we call the two eigenvalues A_ and A+, our new 
problem is to maximize 

Mi,2 = J2\Xk\ = -tX- + {2N-t)X+, (3.39) 
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subject to the constraints 



tA_ + (2iV-t)A+ = 1, (3.40) 

+ a 

~2W 



t\l + {2N -t)\l = ^^^ . (3.41) 



We can now do the problem by solving the constraints for A_ and A+ in terms of t, 
plugging these results into M.i^2-, and then maximizing over t. 

Before continuing, we note two things. First, t cannot be 2N , for if it were, then 
all the eigenvalues would be negative, making it impossible to satisfy Eq. (3.40). 
Second, unless a = 0, t cannot be 0, for if it were, then all the eigenvalues would 
be equal to \/2N by Eq. (3.40), a situation Eq. (3.41) says can occur only if a = 0. 
Since we are not really interested in the case a = 0, for which pa is always the 
maximally mixed state, we assume a > and < t < 2N in what follows. 

Solving Eqs. (3.40) and (3.41) and plugging the solutions into Eq. (3.39), we get 

the two solutions 

N -t± aJt(2N - t) 
-Mi,2 = f^ . (3.42) 

We choose the positive branch, since it contains the maximum. Maximizing with 
respect to t treated as a continuous variable, we obtain the upper bound, 

Mi,2 = Vl + a2 °=^ V2 ~ 1.414 , (3.43) 

which occurs when the degeneracy parameter is given by 

t = N(l , ^ I "g^ 0.292 AT. (3.44) 

The numbers on the right are for the case a = 1, corresponding to the special qubit 
starting in a pure state. Notice that the upper bound (3.43) allows a negativity 
greater than 1 for all a except a = 0. 

Since we did not yet enforce the condition that t be a positive integer, the 
bound (3.43) can be made tighter for specific values of A^ and a by calculating t 
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and checking which of the two nearest integers yields a larger M.i^2- Asymptotically, 
however, the ratio t/N can approach any real number, so this bound for continuous 
t is the same as the bound for integer t in the limit N —>■ oo. 



3.4.2 The s = 1,2,3 Bound 

To deal with this case, we again make the assumption that t of the eigenvalues are 
negative and 2N — t are nonnegative and thus write 

i 2N 

-Mi,2,3 = -$^A,+ Y, A,, (3.45) 

j=i i=t+i 

as before. In addition to the constraints (3.35), we now have a third constraint 

2^ 1 , Q 2 

^l = -^W^ ■ (3.46) 

fc=i 

We specialize to the case a = 1 for the remainder of this subsection, because it is 
our main interest, and the algebra for the general case becomes difficult. 

Introducing three Lagrange multipliers, we can write the function we want to 
maximize as 



t 2N / 2N \ / 2N \ / 2N 

i=l j=t+l \k=l / \k=l / \fc=l 



Ar2 
(3.47) 



Differentiating with respect to Aj and then A^ gives 



-l + fi + 2uXi + 3^X^ = , (3.48) 

+1 + /i + 2iy\j + S^X'^j = . (3.49) 

These equations being quadratic, we see that there are at most two distinct negative 
eigenvalues and at most two distinct nonnegative eigenvalues. Since the sum of the 
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two solutions of either of these equations is — 2z//3.^, however, we can immediately 
conclude either that one of the potentially nonnegative solutions is negative or that 
one of the potentially negative solutions is positive. Hence, we find that at least one 
of the four putative eigenvalues has the wrong sign, implying that there are at most 
three distinct eigenvalues, though we don't know whether one or two of them are 
negative. 

Labelling the three eigenvalues by A, B, and C, we can reduce the problem to 
solving the three constraint equations, 

uA + vB + wC = 1 , 
uA^ + vB^ + wC^ = 1/N , (3.50) 

uA^ + vB^ + wC^ = 1/N^ , 

for A, B, and C and then maximizing A^i,2,3 over the degeneracy parameters u, v, 
and w, which are nonnegative positive integers satisfying the further constraint 

u + v + w = 2N . (3.51) 

We do not associate any particular sign with A, B, and C; the signs are determined 
by the solution of the equations. 

One might hope that the symmetry of Eqs. (3.50) would allow for a simple analytic 
solution, but this appears not to be the case. In solving the three equations, one is 
inevitably led to a sixth-order polynomial in one of the variables, with the coefficients 
given as functions of u, v, and w. Rather than try to solve this equation, which 
appears intractable, we elected to do a brute force optimization for any given value 
of 2N by solving Eqs. (3.50) for each possible value of u, v, and w. Picking the 
solution that has the largest negativity then yields the global maximum. We did this 
for each N up to 2N = 78. The values of m, v, and w that maximize the negativity 
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are always 



u 



--;^) 



1 , w = 2N-l-u, (3.52) 



where [x] denotes the integer nearest to x. The unique eigenvalue corresponding to 
f = 1 is the largest positive eigenvalue, w is the degeneracy of another positive eigen- 
value, and u is the degeneracy of the negative eigenvalue. Notice that the degeneracy 
of the negative eigenvalue is exactly what was found in the s = 1,2 case. Using the 
results (3.52) as a guide, we did a further numerical calculation of the maximum 
for larger values of A^, by considering only the area around the degeneracy values 
given by Eq. (3.52). While this is not a certifiable global maximum, the perturbation 
expansion described below matches so well that the two are indistinguishable if they 
are plotted together for n + 1 > 7. This gives us confidence that the numerically de- 
termined upper bound A^ 1,2,3, which we plot in Fig. 3.7, is indeed a global maximum 
for all N. 

We have used the numerical work to help formulate a perturbation expansion that 
gives the first correction to the N —>■ 00 behavior of the s = 1,2,3 bound. Defining 
X = l/N, we rewrite the constraint equations (3.50) as 

aA + bB + cC = X, 
aA^ + bB^ + cC^ = x\ (3.53) 

aA^ + bB^ + cC^ = x^ 

where a = u/N, b = v/N, and c = w/N. We also have the constraint 

a + b + c = 2 . (3.54) 

As X is the variable that is asymptotically small, we seek an expansion in terms of 
it. 

Our numerical work tells us that there are two positive eigenvalues, one of which 
is larger and nondegenerate. In formulating our perturbation expansion, we let B 
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and C be the positive eigenvalues, with B being the larger one, having degeneracy 
V = bi > 1. We do not assume that 6i is 1, as the numerics show, but rather let 
the equations force us to that conclusion. With this assumption, the form of the 
constraints (3.53) shows that the variables have the following expansions to first 
order beyond the N ^ oo form: 

a = ao + aix^'^ , 

b = bix , (3.55) 

C = Co + CiX^'^ , 

and 

A = Aqx + Aix^^^, 

B = Bix^'^ , (3.56) 

C = Cox + Cix^l^ . 

We see that we are actually expanding in the quantity y = x^^^. In terms of these 
variables, the negativity is given by 

-aA + bB + cC 

Ml,2,3 = 

X 

= -aoAo + cqCo + i-aoAi - aiAo + cqCi + ciCo)x^/^ + 0{x^^'^), (3.57) 

which we now endeavor to maximize. 

Substituting Eqs. (3.55) and (3.56) into the constraints (3.53) and (3.54) and 
equating terms with equal exponents of x, we obtain, to zero order, 

ao + Co = 2 , 
aoAo + coCo = 1 , (3.58) 

(IqAq + CoCq = 1 . 
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Solving for ao, cq, and Cq in terms of Aq and substituting the results into the zero- 
order piece of A^ 1,2,3 gives 



Maximizing Eq. (3.59) gives Aq = 1/2 and, hence, Aq = — 1/-\/2, since A is the 
negative eigenvalue. This leads to oq = 1 — 1/v^, cq = 1 + l/V^, and Cq = 1/v^, 
and the resulting A^ — > 00 upper bound is 7V1 1,2,3 = V^, as expected. 

If we carry this process out to first order beyond the N ^ 00 behavior, we obtain, 
after some algebraic manipulation, A^ 1,2,3 = V2—b^ x^^'^/2'^^^+0{x'^^^). Maximizing 
this simply means making 61 as small as possible, i.e., choosing 61 = 1, whence we 
obtain the following asymptotic expression for the s = 1,2,3 upper bound: 

This shows that the upper bound of v2 is approached monotonically from below 
in the asymptotic regime. In addition, the procedure verifies that in the maximum 
solution, the largest positive eigenvalue is nondegenerate. For the case of qubits we 
have N = 2^, implying that the approach to the N -^ 00 bound is exponentially 
fast. 



3.5 Conclusion 

The mixed-state quantum circuit in Fig. (3.2) provides an efficient method for es- 
timating the normalized trace of a unitary operator, a task that is thought to be 
exponentially hard on a classical computer. If one believes that global entanglement 
is the essential resource for the exponential speedup achieved by quantum computa- 
tion, then the question begging to be answered is whether there is any entanglement 
in the circuit's output state in Eq. (3.3). The purpose of this chapter was to inves- 
tigate this question. 
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A notable feature of the circuit in Fig. (3.2) is that it provides an efficient method 
for estimating the normahzed trace no matter how small the initial polarization a of 
the special qubit in the zeroth register, as long as that polarization is not zero. Since 
all the other qubits are initially completely unpolarized, we are led to characterize 
the computational power of this circuit as the "power of even the tiniest fraction of 
a qubit." We provide preliminary results regarding the entanglement that can be 
achieved for a < 1. Our results are consistent with, but certainly do not demonstrate 
the conclusion that separable states cannot provide an exponential speedup and that 
entanglement is possible no matter how small a is. The question of entanglement 
for subunity polarization of the special qubit deserves further investigation. 

Our key conclusions concern the case where the special qubit is initially pure 
{a = 1). We find that the circuit in Fig. (3.1) typically does produce global entan- 
glement, but the amount of this entanglement is quite small. Using multiplicative 
negativity to measure the amount of entanglement, we show that as the number of 
qubits becomes large, the multiplicative negativity in the state in Eq. (3.1) is a van- 
ishingly small fraction of the maximum possible multiplicative negativity for roughly 
equal splittings of the qubits. This hints that the key to computational speedup 
might be the global character of the entanglement, rather than the amount of the 
entanglement. In the spirit of the pioneering contribution of Wyler [s>>], what hap- 
pier motto can we find for this state of affairs than Multum ex Parvo, or A Lot out 
of A Little. 
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Chapter 4 

Classical simulation of quantum 
computation 



I looked for an answer to my question. But reason could not give me an answer- 
reason is incommensurable with the question. 

-Leo Tolstoy in Anna Karenina 



Progress in our understanding of what makes quantum evolutions computation- 
ally more powerful than a classical computer has been scarce. A step forward, how- 
ever, was achieved by identifying entanglement as a necessary resource for quan- 
tum computational speed-ups. Indeed, a speed-up is only possible if in a quantum 
computation, entanglement spreads over an adequately large number of qubits [52]. 
In addition, the amount of entanglement, as measured by the Schmidt rank of a 
certain set of bipartitions of the system, needs to grow sufficiently with the size 
of the computation [ ]. Whenever either of these two conditions is not met, the 
quantum evolution can be efficiently simulated on a classical computer. These con- 
ditions (which are particular examples of subsequent, stronger classical simulation 
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results based on tree tensor networks (TTN) [75, 20]) are only necessary, and thus 
not sufficient, so that the presence of large amounts of entanglement spreading over 
many qubits does not guarantee a computational speed-up, as exemplified by the 
Gottesman-Knill theorem 



The above results refer exclusively to quantum computations with pure states. 
As shown in Chapter 3, the scenario for mixed-state quantum computation is rather 
different. The intriguing deterministic quantum computation with one quantum bit 
(DQCl or 'the power of one qubit') [ ] involves a highly mixed state that does 
not contain much entanglement [ : ] and yet it performs a task, the computation 
with fixed accuracy of the normalized trace of a unitary matrix, exponentially faster 
than any known classical algorithm. This also provides an exponential speedup 
over the best known classical algorithm for simulations of some quantum processes 
[66]. Thus, in the case of a mixed-state quantum computation, a large amount of 
entanglement does not seem to be necessary to obtain a speed-up with respect to 
classical computers. 

A simple, unified explanation for the pure-state and mixed-state scenarios is 
possible [79] by noticing that the decisive ingredient in both cases is the presence 
of correlations. Indeed, let us consider the Schmidt decomposition of a vector |\I^), 
given by 

X 

1^) = J]Ai|i^)®|iB), (4.1) 

j=i 

where («a| Ja) = (^_b| js) = <^jj and x is the rank of the reduced density matrices 

p^ = trB[|'?/')('?/'|] and Pb = tr^[|^)(^|]; and the (operator) Schmidt decomposition of 

a density matrix p given by [85] 



P 



5^AfO,A®0,s, (4.2) 



i=l 



where tr(Oj^Oj^) = tr(Oj^OjB) = 5ij. The Schmidt ranks x and x^ are a measure 
of correlations between parts A and B, with %" = X^ if p = |\l/)(\l/|. Let the density 
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matrix pt denote the evolving state of the quantum computer during a computation. 
Notice that pt can represent both pure and mixed states. Then, as shown in Refs. 
[79] and [75, 20], the quantum computation can be efficiently simulated on a classical 
computer using a TTN decomposition if the Schmidt rank x^ of p according to a 
certain set of bipartitions A : B oi the qubits scales polynomially with the size of the 
computation. In other words, a necessary condition for a computational speed-up is 
that correlations, as measured by the Schmidt rank xK grow super-polynomially in 
the number of qubits. In the case of pure states (where x = V X^) these correlations 
are entirely due to entanglement, while for mixed states they may be quantum or 
classical. 

Our endeavor in this chapter is to study the DQCl model of quantum computa- 
tion following the above line of thought. In particular, we elucidate whether DQCl 
can be efficiently simulated with any classical algorithm, such as those in [79, 75, 20] 
(and, implicitly, in [52]), which exploit limits on the amount of correlations, in the 
sense of a small x" according to certain bipartitions of the qubits. We will argue here 
that the state pt of a quantum computer implementing the DQCl model displays an 
exponentially large x^ in spite of its containing only a small amount of entanglement 
[17]. We will conclude, therefore, that none of the simulation techniques mentioned 
above can be used to efficiently simulate 'the power of one qubit'. 

On the one hand, our result indicates that entanglement is not behind the (sus- 
pected) computational speed-up of DQCl. On the other hand, by showing the failure 
of a whole class of classical algorithms to efficiently simulate this mixed-state quan- 
tum computation, we reinforce the conjecture that DQCl leads indeed to an expo- 
nential speed-up. We note, however, that our result does not rule out the possibility 
that this circuit could be simulated efficiently using some other classical algorithm. 

Before we move on to the presentation of our results and conclusions, we will 
present a short review of tree tensor networks(TTN) for the sake of completeness 
and continuity. This will be the subject of the next section. 
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4.1 Tree Tensor Networks (TTN) 

The discussion in this section is not original and a more detailed compilation can be 
found in [""]. For simplicity, let us deal with pure states of n qudits, denoted that 
by 

d d 

This state is characterized by c?" complex amplitudes Cj^...j„. We can think of all 
these coefficients as one entity, a rank-ra tensor. Such a tensor can be represented 
graphically, with vertices labelled by c, and connected by n open wires, each of which 
is labelled by a distinct index. One may represent a tensor network by starting with 
such graphical representation of its tensors, and then connecting wires corresponding 
to the same index. As quantum gates are performed on the state |\1/), the graph (or 
the network) gets more involved. 

In the quantum circuit model, result of the quantum computation is obtained by 
taking trace of the final state against measurement operators. In our present picture, 
this translates to contracting vertices of the tensor networks, that is, removing the 
edges between two vertices and replacing them by a single one. This operation is 
known as tensor contraction. The complexity of simulating a quantum computation 
thus reduces to that of simulating tensor contractions. For graphs which have small 
treewdiths, this operation can be approximated classically with a high accuracy in 
polynomial time. Without going into its formal definition, the treewidth of a graph 
tells us how different it is from a tree. A tree has treewidth of 1. Single cycles of 
length at least 3 have a treewidth of 2. 

For any graph G = {V, E), the fine graph G* is defined as follows: V{G*) = E{G) 
and 

E{G*) = {{ei, 62} C E{G) : ei 7^ 62, 3 f G V{G) such that ei and 62 are incident on v} 
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The contraction complexity of a graph G, cc{G) is equal to the the treewidth of 
the line graph G*, tw(G*) [ ]. Although determining the treewidth of a general 
graph is NP-hard [ ], for trees it is trivial. For this simple case then, cc(G') is trivial. 
Additionally, if the graph underlying the quantum circuit is such that tw(G') = d, 
then the circuit can be simulated in time e^^'^\ 

This is the inspiration for defining quantum circuits as tree tensor networks. 
The closer the underlying line graph for a quantum circuit is to a tree, easier it 
is to simulate the corresponding quantum circuit efficiently classically. This, not 
surprisingly then, is the motivation behind all the present simulation techniques of 
quantum computation like matrix product states (MPS), projected entangled-pair 
states (PEPS), Affleck-Kennedy-Leib-Tasaki (AKLT) states. If one is able to show 
that a certain quantum computation cannot be represented as a TTN, then a large 
class of techniques are excluded in providing an efficient classical simulation for it. 
This will be our line of approach in the rest of this chapter, with reference to the 
DQCl model. 



4.2 DQCl and Tree Tensor Networks (TTN) 

The DQCl model, represented in Eq. (3.1), provides an estimate of the normalized 
trace tr(f/„)/2" of a n-qubit unitary matrix f/„ G U(2") with fixed accuracy efficiently 
[55]. This quantum circuit transforms the highly-mixed initial state po = |0)(0| ® 
/„/2'^ at time t = into the final state px at time t = T, 

through a series of intermediate states pt, t G [0,T]. The simulation algorithms 
relevant in the present discussion [52, 79, 75, 20] require that pt be efficiently rep- 
resented with a TTN [75, 20] (or a more restrictive structure, such as a product of 
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fc-qubit states for fixed k [52] or a matrix product state [ ]) at all times t G [0,T]. 
Here we will show that the final state pt, henceforth denoted simply by p, cannot be 
efficiently represented with a TTN. This already implies that none of the algorithms 
in [52, 79, 75, 20] can be used to efficiently simulate the DQCl model. 

Storing and manipulating a TTN requires computational space and time that 
grows linearly in the number of qubits n and as a small power of its rank q. The 
rank g of a TTN is the maximum Schmidt rank xl over all bipartitions Ai : Bi 
of the qubits according to a given tree graph whose leaves are the qubits of our 
system (see [75, 20] for details). The key observation of this chapter is that for a 
typical unitary matrix f/„, the density matrix p in Eq. (4.4) is such that any TTN 
decomposition has exponentially large rank q. By typical, here we mean a unitary 
matrix f/„ efficiently generated through a (random) quantum circuit. That is, f/„ is 
the product of poly(?2) one-qubit and two-qubit gates. In the next section we present 
numerical results that unambiguously suggest that, indeed, typical Un necessarily 
lead to TTN with exponentially large rank q. 

We notice that the results of the next section do not exclude the possibility that 
the quantum computation in the DQCl model can be efficiently simulated with a 
TTN for particular choices of f/„. For instance, if f/„ factorizes into single-qubit 
gates, then p can be seen to be efficiently represented with a TTN of rank 3, and 
we can not rule out an efficient simulation of the power of one qubit for that case. 
Of course, this is to be expected, given that the trace of such Un can be computed 
efficiently in the first place. 



4.3 Exponential growth of Schmidt ranks 

In this section we study the rank q of any TTN for the final state p of the DQCl 
circuit, Eq. (4.4). We numerically determine that a lower bound to such a rank 
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grows exponentially with the number of qubits n. 
The Schmidt rank x oi a. pure state Ip^^^pg) is 

Wa^b) = P \<Pa) \^b) = Y, ^1 <^^^ l*^^) ® <^*s \^b) , (4.5) 

obtained by applying the density matrix p onto a product state |0^) I'j/'b) is a lower 
bound on the operator Schmidt rank x" of p, i.e., x" ^ X- ^^i the purpose of our 
numerics, we consider the pure state f/„ |0) ". We build t/„ as a sequence of 2n 
random two-qubit gates, applied to pairs of qubits, also chosen at random. The 
random two-qubit unitaries are generated using the mixing algorithm presented in 
[30]. Note that applying 2n gates means that the resulting unitary is efficiently 
implementable, a situation for which the DQCl model is valid. For an even number 
of qubits n, we calculate the smallest Schmidt rank x over all n/2 : n/2 partitions 
of the qubits (similar results can be obtained for odd n). The resulting numbers are 
plotted in Fig (4.1). 

The above numerical results strongly suggest that the final state p in the DQCl 
circuit has exponential Schmidt rank for a typical unitary f/„. We are not able to 
provide a formal proof of this fact. This is due to a general difficulty in describ- 
ing properties of the set I[Jgc(2") of unitary matrices that can be efficiently realized 
through a quantum computation. Instead, the discussion is much simpler for the set 
I[J(2") of generic n-qubit unitary matrices, where it is possible to prove that p cannot 
be efficiently represented with a TTN for a Haar generated f/„ G U(2"), as discussed 
in the next section. Notice that Ref. [28] shows that random (but efficient) quantum 
circuits generate random ra-qubit gates f/„ G lJqc{2"') according to a measure that 
converges to the Haar measure in U(2"). Combined with the theorem in the next 
section, this would constitute a formal proof of the otherwise numerically evident 
exponential growth of the rank q of any TTN for the DQCl final state p. 



Chapter 4. Classical simulation of quantum computation 



82 



Lower bound on the Schmidt rank of the DQC1 state 




Number of qubits 

Figure 4.1: Lower bound for the operator Schmidt rank x'^ of the DQCl state for 
any equipartition n/2 : n/2, as given by the Schmidt rank x of the pure state in 
Eq. (4.5). The dots are for even numbers of qubits, and the fit is the hne 2"/^. x 
is calculated for a pure state obtained by applying 2n random 2-qubit gates on the 
state |0) ^. This is evidence that for a typical unitary Un, the rank q of any TTN 
for the DQCl state p in Eq. (4.4) grows exponentially with n. 



4.4 A formal proof for the Haar-distributed case 

Our objective in this section is to analyze the Schmidt rank x" of the density matrix 
p in Eq. (4.4) for certain bipartitions of the n + 1 qubits, assuming that Un G U(2"') 
is Haar-distributed. As an aside, we calculate the negativity of a random pure state 
in Appendix C. 



It is not difficult to deduce that for any tree of the n + 1 qubits, there exists at 
least one edge that splits the tree in two parts A and B, with n^ and ub qubits, 
where Uq = min(n^,nB) fulfills n/5 < tiq < 2n/5. In other words, if a rank-g TTN 
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exists for the p in Eq. (4.4), then there is a bipartition of the n + 1 qubits with uq 
qubits on either Aoi B and such that the Schmidt rank x^ ^ Q- Theorem 6, our main 
technical result, shows that if f/„ is chosen randomly according to the Haar measure, 
then the Schmidt rank of any such bipartition fulfills x" ^ 0(2"°). Therefore for a 
randomly generated f/„ G 0(2"'), a TTN for p has rank q (and computational cost) 
exponential in n, and none of the techniques of [52, 79, 75, 20] can simulate the 
outcome of the DQCl model efficiently. 

Consider now any bipartition A : B oi the n + 1 qubits, where A and B contain 
ua and ub qubits, with the minimum uq of those restricted by n/5 < no < 2n/5. 
Without loss of generality we can assume that the top qubit lies in A. Actually, we 
can also assume that A contains the top ua qubits. Indeed, suppose A does not have 
the Ua top qubits. Then we can use a permutation P„ on all the n qubits to bring 
the Ua qubits of A to the top ua positions. This will certainly modify p, but since 

P„ J \u„ I„ J \ P^ J \ v„ /„ j 

where Vn = PnUnPn is another Haar-distributed unitary, we obtain that the new 
density matrix is of the same form as p. Finally, in order to ease the notation, we 
will assume that ua = n^ (identical results can be derived for ub = no). Thus 
n/5 < n^ < 2n/5. 

We note that 
L Ul 



'n 




l2®In+ I \ ®Ul+ \ I ®f/„, (4.7) 

SO that if we multiply p by the product state 

\<t)s) = \t,t,j) = \t,tA)\3B), (4.8) 

where a = {t,i,j), t = 0,1; i = 1, .. . cIa] j = 1, . . . cIb, we obtain \ip^) = p 10^) where 

Wd) = { ^ ^ (4.9) 

2^(|l,z,j) + |0)®?7t|z,j)) ift = l. 
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This also justifies our choice of the pure state used in the numerical calculations in 
the previous section. 

Let us consider now the reduced density matrix 



R 



= tr A[\i^d){i^6 



= ^A\3)(:j\+^'A[Un\^,J){^,J\Ul]) (4.10) 

for t = (for t = 1, Un and U^ need to be exchanged). For a unitary matrix f/„ 
randomly chosen according to the Haar measure on U(n), f/„ \i,j) is a random pure 
state on A® B. Here, and henceforth A is the space of the first ua qubits without 
the top qubit. It follows from [43] that the operator 

Q = tTA[Un\l,j){t,j\Ui] (4.11) 

has rank (Ia- Therefore the rank of a? (equivalently, the Schmidt rank x of iV^a)) 
is at least 2"°. From Eq. (4.5) we conclude that the Schmidt rank of p fulfills 
x" > 2"" > 2"/^. We can now collate these results into 

Theorem 6. Let Un be an n-quhit unitary transformation chosen randomly according 
to the Haar measure on U(2"), and let A : B denote a bipartition ofn + 1 qubits into 
Ha andriB qubits, where Uq = min(nA,'n-B)- Then n/5 <nQ< 2n/5 and the Schmidt 
decomposition of p in Eq. (4-4) according to bipartition A : B fulfills x^ > 2"/^. 

We have seen that we cannot efficiently simulate DQCl with an algorithm that 
relies on having a TTN for p with low rank q. However, in order to make this 
result robust, we need to also show that p cannot be well approximated by another 
p accepting an efficient TTN. This is what we do next. 

We now explore the robustness of the statement of Theorem 6. To this end, we 
consider the Schmidt rank x" for a density matrix p that approximates p according 
to a fidehty F{Oi, O2) defined in terms of the natural inner product on the space of 
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linear operators, 



F(0i,02) = tr(0l02) / Vtr(0l0i)^tr(0|02) , 

where F = 1 if and only if Oi = O2 and F = \ {ipil 1P2) P for projectors Oj = P^. on 
pure states \ipi). We will show that if p is close to p, then x^ for a bipartition as in 
Theorem 1 is also exponential. To prove this, we will require a few lemmas which 
we now present. 

Lemma 2. Let |\E') he a bipartite vector with x terms in its Schmidt decomposition, 

X X 

where N,^ = ^y{W\W), and let |$) be a bipartite vector with norm N^ and Schmidt 
rank x' , where x' ^ X- Then, 



max I (M/l $) I = N^N^ 



\ 



EA2. (4.12) 



1=1 



Proof: Let /ij denote the Schmidt coefficients of |$). It follows from Lemma 
1 in ["'] that max|$) | (^|$) | = Ni!,Nq,^f^^ Xifii, and the maximization over /ij 
is done next. A straightforward application of the method of Lagrange multipliers 
provides us with /Xj = cAj, z = 1, 2, . . . , x' for some constant c. Since ^2^=1 ^A — ^ — 



c^Ef=iA?,c=l/^Ef=iA?.Thus, 

x' 
max I (^1 $) I = cN^N^ V A? 

1=1 

and the result follows. D 

We will also use two basic results related to majorization theory. Recall that, by 
definition, a decreasingly ordered probabihty distribution p = {pi^p2, . . . ,Pd), where 
Pa > Pa+i > 0, X^aPa = 1, IS majouzcd by another such probability distribution q, 
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denoted p ~< q, ii q is more ordered or concentrated than p (equivalently, p is flatter 
or more mixed than q) in the sense that the following inequalities are fulfilled: 

k k 

5^Pa<5^ga yk=l,...,d (4.13) 

with equality for k = d. The following result can be found in Exercise II. 1.15 of [7]: 

Lemma 3. Let pg and p^ he density matrices with eigenvalues given by probabil- 
ity distributions x and y. Let cr(M) denote the decreasingly ordered eigenvalues of 
hermitian operator M. Then 

cr{px + Py) ^x + y. 

The next result follows by direct inspection. 

Lemma 4. Let coefficients Si, 1 < i < d, be such that —5 < 6i < S for some positive 
6 <1 and Yli^i — 1; ^'^^ consider the probability distribution p{{6i}) , 
^.rx^^ fl l + Si I + 62 I + 6d 



2 2rf ' 2rf ' ' 2d 
Then 

v{{5^})-<pm}). 

where 

5 i<d/2 



-6 i > d/2, 



and we assume d to be even. 



Finally, we need a result from [43]: 
Lemma 5. With probability very close to 1, 



Pr 



:i-5)-<Q<il + S)- 
ClA dA 



> 1 _ ( 12_^^ "" 2^~'^B <57l41n2) 

>^-0{ ../ I .^ , (4.14) 

exp(o"^exp(n)j/ 
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where cLa = 2"-* = 2"° and cLb = 2"'-s = 2"""""'"^, and the operator Q defined in Eq. 
(4-11) is within a ball of radius 6 of a (unnormalized) projector T/dA of rank dA 
[provided ds is a large multiple o/ d^logd^/^^ [ ], which is satisfied for large n, 
given that n/5 < Uq < 2n/5]. 

Our second theorem uses the fact that the Schmidt decomposition of p does not 
only have exponentially many coefficients, but that these are roughly of the same 
size. 

Theorem 7. Let p, Un, and A:B be defined as in Theorem 6. IfF{p,p) > 1 — e, then 
with probability p{S,n) = 1 — 0(exp(— (5^exp(n))), the Schmidt rank for p according 
to bipartition A:B satisfies x" > (1 — 4e — (5)2"/^. 

Proof: For any product vector of Eq. (4.8) we have 



{ttj\pp\ttj)\ < N^N, 

< N^N^gix^/dA 



vlt 



X' 



X E(^l')' (4.15) 

\ k=i 



where 



, 1 + (1 + 6)x , , 

9i^)^\l 4~^' ^^-^^^ 

and N^ = \/ {tij\ p^ \tij), iV^ = ^Jitijl p^ \tij). The first inequality in (4.15) follows 
from Lemma 1, whereas the second one follows from the fact that the spectrum p of 

PB = {N^YhiA[p\timij\p\ = ^(|j)(j| + g), 

where Q has all its dA non-zero eigenvalues g^ in the interval 2^"o(l — 5) < Qi < 
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2~'^°(1 + 6), is majorized by p{{S*}), as follows from Lemmas 2 and 3. Then, 

tipp 



l-e< 



A/trp^A/trp^ 

E5(«IPPl") 






< 9ix^/dA), 

where in the last step we have used the Cauchy-Schwarz inequality, | (x||/) | < 
\/{x\ x)\J {y\ y). The result of the theorem follows from g{x^ 12"^°) > 1 — e. D 



4.5 Conclusions 

The results in this chapter show that the algorithms of [52, 79, 75, 20] are unable 
to efficiently simulate a DQCl circuit. The efficiency of a quantum simulation using 
these algorithms relies on the possibility of efficiently decomposing the state p of the 
quantum computer using a TTN. We have seen that for the final state of the DQCl 
circuit no efficient TTN exists. 

It is also interesting to note that the numerics and Theorems 1 and 2 in this 
chapter can be generalized for any fixed polarization r (0 < r < 1) of the initial 
state t|0)(0| + (1 — t)1/2 of the top qubit of the circuit in Eq (3.1), implying that 
the algorithms of [52, 79, 75, 20] are also unable to efficiently simulate the power of 
even the tiniest fraction of a qubit. 
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Chapter 5 

Quantum discord in the DQCl 
model 



We are coining now rather into the region of guesswork! 
Say, rather, into the region where we balance probabihties and choose the most 
hkely. It is the scientific use of imagination, but we always have some material basis 
on which to start our speculation. 

-Arthur Conan Doyle in The Hound of the Baskervilles 



We have seen in the last two chapters (Chapters 3, 4) that the DQCl model 
presents to us a very intriguing challenge. It is a mixed-state quantum computation 
scheme [nri] that has very little entanglement [17], yet it cannot be simulated by 
matrix product state techniques [j -^j. This suggests that entanglement cannot be 
the sole resource that drives mixed-state quantum computation. Although in pure 
state quantum computation entanglement can be shown to be an essential resource 
[52, 79], mixed-state quantum computation is a different story. It is evident that there 
is a gap in the resource based accounting for quantum computational speedups. It 
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is this gap that we hope to fill with the quantum discord [!;']. 

Quantum discord, introduced by Ollivier and Zurek [ ], captures the nonclassi- 
cal correlations, including but not limited to entanglement, that can exist between 
parts of a quantum system. Some details about quantum discord are presented in 
Sec 2.2. One way of studying the quantum nature of a computational process is 
to investigate the nonclassical correlations in the quantum state at various stages 
during the computation. We investigate the effectiveness of discord in characteriz- 
ing the performance of the model of quantum information processing introduced by 
Knill and Laflamme in [55], which is often referred to as the power of one qubit, 
or DQCl. In this model, information processing is performed with a collection of 
qubits in the completely mixed state coupled to a single control qubit that has some 
nonzero purity. Such a device can perform efficiently certain computational tasks for 
which there is no known efficient method using classical information processors. 

In this thesis, thus far we have seen how discord can be used to characterize the 
nonclassical nature of the correlations in quantum states and studied its role and 
place in quantum information theory. We now apply these ideas to the DQCl or 
power- of- one- qubit model [' ] of mixed-state quantum computation, which accom- 
plishes the task of evaluating the normalized trace of a unitary matrix efficiently. 
The quantum circuit corresponding to this model has a collection of n qubits in the 
completely mixed state, In/'^"', coupled to a single pure control qubit. A generalized 
version of this quantum circuit, with the control qubit having sub- unity polariza- 
tion is shown below in Fig 3.1. This circuit evaluates the normalized trace of Un, 
T = tT{Un)/2^, with a polynomial overhead going as 1/a^. 

The problem of evaluating r is believed to be hard classically. Quantum mechan- 
ically, the circuit provides an estimate of r up to a constant accuracy in a number 
of trials that does not scale exponentially with n. It does so by making X and Y 
measurements on the top qubit. The averages of the obtained binary values provide 
estimates for tr = Re(r) and r/ = Im(r). The top qubit is completely separable 
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from the bottom mixed qubits at all times. The final state has vanishingly small 
entanglement, as measured by the negativity [ ] across any split that groups the 
top qubit with some of the mixed qubits. Nonetheless, there is evidence that the 
quantum computation performed by this model cannot be simulated efficiently using 
classical computation [18]. 

The DQCl circuit transforms the highly-mixed initial state po = |0)(0| ® /n/2" 
into the final state pn+i given by Eq 3.1. Within this model the only place to look 
for nonclassical correlations is in this state. 

Everything about the DQCl setup, including the measurements on the control 
qubit, suggests a bipartite split between the control qubit M and the mixed qubits 
S. Relative to this split, we turn to computing the quantum discord for the state 
PsM = Pn+i- The joint state pn+i has eigenvalue spectrum 

A(p„+i) = 7^— t(1 - a, ■ ■ ■ , 1 - a, 1 + a, ■ ■ ■ , 1 + a), 

2"times 2"times 

which gives a joint entropy H{S, M) = n + -f^2[(l — ct)/2], where H2[-] is the binary 
Shannon entropy. The marginal density matrix for the control qubit at the end of 
the computation is 



PM = 7^ , (5.1) 




which has eigenvalues (1 ± a;|r|)/2 and entropy H{M) = -f^2[(l ~ Q^kl)/2]- 

The evaluation of the quantum conditional entropy involves a minimization over 
all possible one-qubit projective measurements. The projectors are given by II-i- = 
l(/i ± a • cr), with a • a = a^ + Og + Og = 1. The post-measurement states are 

Ps\± = ^^;^ i^In ± a^ f/n ± «^^— f^IJ ' (5-2) 

occurring with outcome probabilities p-t = [l±a{aiT^+a2Tj)]/2. The post-measurement 
states are independent of as, so without loss of generality, let a^ = 0, ai = cos 0, and 
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02 = sin (p. The corresponding post-measurement states are 

Ps\± = ^^^. \^n ± a ^ J . (5.3) 

To find the discord of the state at the end of the computation, we need the 
spectrum of ps\± so that we can compute H{ps\±). The eigenvalues of any unitary 
operator f/„ are phases of the form 6*^*= , so we have 

A J ^^ ^j =cos(^fc-0), A; = l,---,2", (5.4) 

and 

1 l±acos(^fe-0) 

Afc(p5|±) = 7^ ^ , , — —TV = qk±- (5.5) 

2" 1 ± a(rij cos + r/ sm 0) 

We also have tr = 2~"^^cos6'fc and tj = 2^" ^^sin6'yfc. All this gives H{ps\±) = 
H{q±) and thus 

Hu± = p+H{ps\+)+P-H{ps\^) 

1 a 

= -[H{q+) + H{q^)] + -(r^,cos</) + r, sin0)[ff (g+) - i7(g_)]. (5.6) 

We now use the fact that we are interested in the behavior of the quantum 
discord of the DQCl state for a typical unitary. By typical, we mean a unitary 
chosen randomly according to the (left and right invariant) Haar measure on 0(2"'). 
For such a unitary, it is known that the phases 9k are almost uniformly distributed 
on the unit circle with large probability [ ]. Thus for typical unitaries ^^ 6*^*= is 
close to zero. Hence both r^j and tj are small, and we can ignore the second term on 
the right-hand side in Eq. (5.6). In addition, the phases 6k can be taken to be placed 
at (with large probability) the 2"'th roots of unity, i.e., 9k = 27rfc/2". It follows that 
the spectra Afc(p5'|-|-) are independent of 0. Hence the entropies we are interested 
in computing are also independent of (p, and we can set (p to zero without loss of 
generality. This choice for corresponds to measuring the pure qubit M along X. 
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The X measurement gives the real part of the normahzed trace of Un, and it is one 
of the two measurements discussed in the original proposal by Knill and Laflamme. 
Setting (f) = 7t/2 yields the other measurement, along Y, which gives the imaginary 
part of the normalized trace of Un- 

In the limit of large n, we can simplify Eq. (5.6) as follows: 



H = -[H{q+) + Hiq_)] 



nn 

— T 



fc=l 



1 + Q;cos6'fc)log 



1 + a cos 6k 



2" 



+ (1 — a;cos6'jt) log 



\ — a cos 6*^ 



n 



on 

-—y 



k=\ 



log (l — C? COS^ Qk) + a COS ^fc log 



1 + a COS Q}^ 
1 — a cos 6k 



■ (5.7) 



Furthermore, when n is large, we can replace the sum in the above equation with an 
integral to obtain 



H 



n — 



An 



1 + a cos X 



2iT p2n 

log(l — a^ cos^ x)dx + a cos x log ( — '- 1 dx 

Jo V 1 — a cos X 



n 



+ 1 - log (l + Vl - aA - (l - Vl - aA loge. 



(5.8) 



Note that when the sums are replaced by integrals, H{q^) — H{q_) = 0, providing 
further justification for ignoring the second term in Eq. (5.6). 

When \t\ is small, H{M) ~ 1, and the quantum discord for the DQCl state is 
then given by the simple expression 

'1 — a" 



V 



DQCl 



2 - ^2 (^^) - log (l + VT^^) - (l - v^r^^) log e. (5.9) 



Note that the above expression for the discord, valid for large n. is independent 
of n. Figure 5.1 compares the discord from Eq. (5.9) with the average discord in a 
DQCl circuit having five qubits in the mixed state (n = 5) coupled to a control qubit 
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with purity a. The average is taken over 500 instances of pseudo-random unitary 
matrices, generated using the efficient algorithm presented in [■-,'..]. The convergence 
of this ensemble to the Haar measure on the unitary group is shown in [j.s\. We 
see that in spite of the approximations made in obtaining Eq. (5.9), the analytic 
expression provides a very good estimate of the discord even when n is as low as five. 



Q 0.3 




Figure 5.1: The dashed (red) line shows the average discord in a DQCl circuit 
with five qubits in the mixed state (n = 5) coupled to a qubit with purity a. The 
average is taken over 500 instances of pseudo-random unitary matrices. The discord 
is shown as a function of the purity of the control qubit. The solid (green) line shows 
the analytical expression in Eq. (5.9), which grows monotonically from at a = 
(completely mixed control qubit) to 2 — loge = 0.5573 at a = 1 (pure control qubit). 
These values of discord should be compared with a maximum possible discord of 1 
when M is a single qubit. 



There is no entanglement between the control qubit and the mixed qubits in the 
DQCl circuit at any point in the computation, yet there are nonclassical correlations, 
as measured by the discord, between the two parts at the end of the computation 
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for any a > 0. Other bipartite splittings of p„+i can exhibit entanglement, but it 
was shown in [I ~] that the partial transpose criterion failed to detect entanglement 
in pn+i for tt < 1/2. In this domain, several other tests for entanglement, including 
the first level of the scheme of Doherty et al. [ ], which is based on semi-definite 
programming, also failed to detect entanglement. The above expression is thus the 
first signature of nonclassical correlations in the DQCl circuit for a < 1/2. 

In conclusion, we calculated the discord in the DQCl circuit and showed that 
nonclassical correlations are present in the state at the end of the computation even 
if there is no detectable entanglement. Thus for some purposes, quantum discord 
might be a better figure of merit for characterizing the quantum resources available to 
a quantum information processor. We present evidence of the presence of nonclassical 
correlations in the DQCl circuit when a < 1/2. For qubits quantum discord is known 
to be a true measure of nonclassical correlations [11]. This suggests that nonclassical 
correlations other than entanglement, as quantified by the discord, might explain the 
(sometimes exponential) speed-up in the DQCl circuit and perhaps the speedup in 
other quantum computational circuits. For pure states, discord becomes a measure 
of entanglement. Therefore, using discord to connect quantum resources to the 
advantages offered by quantum information processors has the additional advantage 
that it works well for both pure- and mixed-state quantum computation. 
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Chapter 6 



Conclusion 



If you want a happy ending, that depends, of course, on where you stop your story. 

-Orson Welles 



The intent of this chapter is eponymous. The attempt will be to present the 
results of this thesis in a well-rounded manner, put them all in the perspective of a 
'big picture', and end with a few suggestions for open problems. 

A one line conclusion of this thesis would be the following: 

There is more to quantum information science than just entanglement. 

However, this statement needs qualifications. I presented evidence that entangle- 
ment fails to explain exponential speedups in mixed-state quantum computation, at 
least in the DQCl model of Knill and Lafiamme, the reason being that the amount 
of entanglement in the system is minimal, not scaling with the system size. This was 
the content of Chapter 3. Then I showed that this lack of entanglement does not 
mean that the system is classically simulatable. This is important as the simulation 
techniques that have their roots in renormalization group theory are very powerful 
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and among the most versatile in many-body physics. Thus Chapter 4 leaves us with 
a deep foreboding about the role of entanglement in mixed-state quantum compu- 
tation and the criterion for classical simulatability of quantum systems that involve 
entanglement. In Chapter 5, I was able to show that there is a non-trivial amount 
of quantum discord in the DQCl circuit. It is thus possible that quantum discord is 
the resource that drives mixed-state quantum computation. 

Quantum discord can be thought of as a generalization of entanglement. For 
pure states, they are identical. For mixed quantum states, discord captures non- 
classical correlations beyond entanglement. As I have presented in Section 2.3, the 
scope of quantum discord is wider than just mixed-state quantum computation. I 
discussed the role of discord in distribution of entanglement without investing any 
entanglement. 

Every PhD thesis is a compilation of results and solutions to problems that often 
seem disparate. Yet, it is rarely the case that these problems are not steps towards 
comprehending a fundamental problem in the field of study. In my case, the results 
dealt with the DQCl model, its entanglement content, its classical simulatablity 
or resource behind its exponentially enhanced operation. I also proposed quantum 
discord as a substitute for entanglement for the role of the resource in mixed-state 
quantum computation. I also studied how discord seems to be relevant in other 
information processing tasks and the major part that the DQCl model can play in 
bringing together different areas of physics, mathematics and computer science. The 
fundamental problem I have attempted to attack is the following: 

Is entanglement a key resource for computational power? 

This is what I would call the 'big picture.' Jozsa and Linden have spoken most lucidly 
on this. The answer to this proceeds as follows: The significance of entanglement 
for pure-state computations is derived from the fact that unentangled pure states of 
n qubits have a description involving only poly{n) parameters (in contrast to 0(2") 
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parameters for a general pure state). But this special property of unentangled states 
(of having a 'small descriptions) is contingent on a particular mathematical descrip- 
tion, as amplitudes in the computational basis. If we were to adopt some other 
choice of mathematical description for quantum states (and their evolution), then, 
although it will be mathematically equivalent to the amplitude description, there will 
be a different class of states which will now have a polynomially sized description; 
i.e. two formulations of a theory which are mathematically equivalent (and hence 
equally logically valid) need not have their corresponding mathematical descriptions 
of elements of the theory being interconvertible by a polynomially bounded compu- 
tation. With this in mind we see that the significance of entanglement as a resource 
for quantum computation is not an intrinsic property of quantum physics itself, but 
is tied to a particular additional (arbitrary) choice of mathematical formalism for 
the theory. 

Thus, suppose that instead of the amplitude description we choose some other 
mathematical description V of quantum states (and gates). Indeed, there is a rich 
variety of possible alternative descriptions. Then there will be an associated property 
of states, prop{V), which guarantees that the D-description of the quantum compu- 
tational process grows only polynomially with the number of qubits (e.g. if V is the 
amplitude description, then prop{T>) is just the notion of entanglement). Thus, just 
as for entanglement, we can equally well claim that prop{T>) for any T) is an essential 
resource for quantum-computational speed-up! Entanglement itself appears to have 
no special status here. 

Explicit examples of prop{T>) would include the dimension of the Hilbert space 
in which the state resides, the purity of a quantum state or the stabilizer formalism 
and even quantum discord. Thus, in a fundamental sense, the power of quantum 
computation over classical computation ought to be derived simultaneously from all 
possible classical mathematical formalisms for representing quantum theory, not any 
single such formalism and associated quality (such as entanglement), i.e. we have 
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arrived at the enigmatic prospect of needing a representation of quantum physics 
that does not single out any particular choice of mathematical formalism. 

We now seem to have realized the fundamental problem we ought to be tackling is 
perhaps too hard to solve. Yet, this should not be a detriment, but rather a challenge 
to be at least working on the right problem, hard or not . In this thesis, my attempt 
has been to start chipping at a small corner of this massive monolith. The sculpture 
is far from being done, but the work has begun. 

As is with any work of scientific research, there always remain fairly immediate 
questions whose answers hold the possibility of further enlightening the topic at 
hand. This thesis is no exception. There remain open questions that could have 
been addressed in the course of my research. It is my belief that given enough time 
and effort, these could be resolved. Some of these concern the DQCl model, others 
the nature and properties of quantum discord. As the last segment of this thesis, I 
sketch some of these prospective problems: 

1. DQCl and knot theory: Evaluation of the Jones polynomial of the trace 
closure of a braid is DQCl complete. Knot invariants are related to partition 
functions of the Ising model, and this can be used to find a physical interpre- 
tation of the complexity class DQCl. As a continuation of the discussion in 
Section 2.2, the actual task would be to cast the complexity of approximat- 
ing partition functions with different boundary conditions as graph theoretic 
problems whose complexity is well studied. 

2. Classical simulation of quantum systems: In Chapter 4, I found that the 
DQCl system is not efficiently simulatable classically using matrix product 
states (MPS) which had successfully simulated other quantum systems with 
little entanglement [77, 85]. It is known that pair entangled projected states 
(PEPS) can simulate any quantum computation efficiently [70]. My next aim 
in this regard is to study the simulatability of the DQCl model using some 
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intermediate simulation scheme, lying between MPS and PEPS, and if pos- 
sible, find a measurement-based quantum computation system that emulates 
DQCl. This will tell us more about the limits of quantum computation and 
the boundaries of classical computation. 

3. DQCl and discrete Wigner functions: Quantum computational speedup 

has been related to the positivity of discrete Wigner functions [ ' ']. Such 
Wigner functions W can be defined so that the only pure states having non- 
negative W for all such functions are stabilizer states. Also, the unitaries pre- 
serving non-negativity of W for all definitions of W in the class form a subgroup 
of the Clifford group. This means pure states with non-negative W and their 
associated unitary dynamics are classical in the sense of admitting an efficient 
classical simulation scheme using the stabilizer formalism [11]. It is not hard to 
calculate discrete Wigner functions for mixed states. The hindrance in draw- 
ing conclusions similar to the pure state case arises from the non-uniqueness 
of discrete Wigner functions as they now are defined [•)•'»]. One must mini- 
mize over all such definitions to draw conclusions about the non-negativity of 
Wigner functions as is done in [II]. Further thought in this direction might be 
worthwhile and provide new insights into the power of the DQCl model. 

4. Computation of quantum discord: The minimization involved in the com- 
putation of quantum discord is non-trivial. It is a minimization of an entropic 
quantity over the set of POVMS. This quantity can be recast as the Holevo 
quantity in quantum information theory. We know that the computation of 
Holevo capacity in general is NP-Hard [ ], but the quantum conditional en- 
tropy is a very special case. An interesting question would therefore be the 
complexity of calculating the quantum conditional entropy. More practically, 
attempts to reduce the minimization to a convex optimization technique, like 
semi-definite or semi-infinite programming, would be helpful and worthwhile. 

5. Quantum discord in quantum communication: Beyond the results of 
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Section 2.3, the most important contribution would be to prove the role of 
discord in the general protocol of entanglement distribution. The protocol is 
given in [15], but the challenge is that the general protocol involves a mediating 
particle which is higher than two dimensional, making the explicit evaluation 
of discord challenging. In this context, the previous task is related to this one. 

The above list is by no means exhaustive. The endeavor has been to present a few 
problems that I myself have been interested in. These might lead us to understanding 
certain aspects of the DQCl model and quantum discord. It is likely that their 
study would lead to more open problems. All these and many more challenging and 
profound questions can arise in the study of mixed-state quantum computation and 
quantum discord. 
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Appendix A 



DQCl and Knot Theory 



A.l Quadratically Signed Weight Enumerators 

A general quadratically signed weight enumerator is of the form 

S{A,B,x,y)= J2 (-l)'"'''a:IS"-l'l, 

b:Ab=0 

where A and B are be matrices over Z2 with B of dimension n hj n and A of 
dimension m by n. The variable b in the summand ranges over column vectors in 
Z2 of dimension n, b^ denotes the transpose of b, \b\ is the weight of b (the number 
of ones in the vector b), and all calculations involving A, B and b are modulo 2. 
The absolute value of S{A,B,x,y) is bounded by {\x\ + \y\)"'. In general, one can 
consider the computational problem of evaluating these sums. In particular, it is 
known that for integers k, I, evaluating S{A, B, k, I) exactly is ^^P complete. Now, 
let A be square, of size n by n, and let A denote the lower triangular part of A, 
which is the matrix obtained from A by setting to zero all the entries on or above 
the diagonal. Let diag{A) denote the diagonal matrix whose diagonal is the same as 
that of A and I denote the identity matrix. For matrices C and D with the same 
number of columns, let [C; D] denote the matrix obtained by placing C above D. 



Appendix A. DQCl and Knot Theory 104 

Then, the following two theorems hold: 

Theorem 8. Given diag{A) =1, k,l positive integers, and the promise that 

\S{A, A,k,l)\ > {k'^+P)"-^^/2, determining the sign ofS{A, A,k,l) is BQP-complete. 

Theorem 9. Given diag{A) =1, k,l positive integers, and the promise that 
\S{[A;A'^],A,k,l)\ > (A;2 + /2)n/2/2^ the sign of S{[A; A^],'A,k,l) can be efficiently 
determined by the DQGl model. 

In these two problems, the integers k and / can be restricted to 4 and 3, respec- 
tively, without affecting their hardness with respect to polynomial reductions (using 
classical deterministic algorithms). Note that in the second one, the question of 
DQCl completeness is still open. 

Lidar has made the first connection between QSWEs and partition functions of 
spin models on certain classes of graphs. In that case, A denotes the adjacency ma- 
trix of the relevant graph. Though these classes are highly restricted at the present, 
attempts at broadening their scope are ongoing. Through the connection to QSWEs, 
DQCl can be applied to analyze problems in coding theory, like weight generating 
functions of binary codes. Studying the complexity of such problems will certainly 
give us a better understanding of the complexity classes P, BQP, and DQCl and 
their relative inclusions. QWSEs relate the DQCl model to Ising model partition 
functions in a way entirely independent from that to be outlined in Section A. 2. On 
the other hand, as shown in Fig 1.2, the DQCl model can be employed to tackle 
problems in graph and knot theory, disciplines that have provided scores of problems 
for theoretical computer science in general and quantum computation in particular. 
Identifying problems that are DQCl complete and studying then in relation to BQP, 
P and NP complete graph and knot theoretic problems might unravel the underly- 
ing mathematical structures that separate P, BQP, DQCl. Such an understanding 
might illuminate greatly the mathematical framework behind the power of quantum 
computation, just as studying the complexity of partition functions might unravel 
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the physical foundations for that power. The diagram in Fig 1.2 outhnes prospective 
avenues for exploring such problems in the context of mixed-state quantum com- 
putation, and has the potential of opening up whole new paradigms in mixed-state 
quantum computation research. 

A. 2 Jones Polynomials 

It has been known for quite some time that evaluating the partition functions of 
certain spin models, like the Ising model, is related to knot invariants, in particular, 
the Jones polynomial [ ']. This is independent of the connections arrived at via 
Kauffman brackets and QSWEs [59, 56] as we discussed in Appendix A. The differ- 
ence between trace and plat closure is in the way the ends of the braids are joined to 
make a knot. This has a counterpart in the physical world. Plat closures correspond 
to open boundary conditions, while trace closures correspond to periodic boundary 
conditions [50]. The partition function of a system depends on its boundary condi- 
tions. Therefore, one would expect that its evaluation would behave analogously. Of 
course, when we talk of evaluation, we actually mean approximation to its true value 
with a reasonable accuracy. The notion of 'reasonable accuracy' is a fairly technical 
one. Put simply, we aim for an approximation whose error scales only polynomially 
with the problem size. The exact evaluation of partition function of the Ising model 
on a general graph is T^P-hard even when all the interactions are unity and there 
are no external magnetic fields [81]. We believe that the exact evaluation of the 
normalized trace of a unitary is #P-complete (See Section 3.1). As the DQCl model 
provides a polynomial algorithm for an approximate solution to this hard problem, 
it might be possible to approximate partition functions through the DQCl algorithm 
and uncover physical reasons for the success of the power of one qubit model. 
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Appendix B 



Proof of the Lemma 



Lemma: tr(Afi) = ti^AB). 



Proof: Define two operators A and B by 

A='^aij^ki\ij){kl\ , ^ = X^ bmn,pq\'mn){pq\ . (B.l) 

i,j,k,l m,n,p,q 

Taking the partial transpose with respect to the second subsystem, we find 

A='^aij^kiW){kj\ , ^ = X^ bran,pq\'mq){pn\ . (B.2) 

i,j,k,l m,n,p,q 

Calculating the quantities of interest, we find that they are indeed equal. 

tr{AB) = ^ aij^kibmn,pq{pn\il) {kj\mq) = '^ aij^kibki,ij , (B.3) 

m, n,p, q 

tr (AB) = ^ aij^kibmn,pq{pq\ij) {kl\mn) = '^ aij^kibki,ij ■ (B.4) 

i,j,k,l, '■" ' 

m,n,p, q 
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Appendix C 



Negativity of a random pure state 



An eigenvector of a random n-qubit unitary can be considered to be a random pure 
state in the Hilbert space of n-qubits, Ti. Any such state can be Schmidt decomposed 
for a given bipartition 1-La®'Hb with dimensions /i and v {^ <v) and the distribution 
of the Schmidt coefficients is given by (for n = u) 

j=l l<i<j<tJ. k=l 

Since the negativity for pure states is 



r ^ 



(C.3) 






At this point, it helps to change variables such that Pi = rqi which removes the 
hurdle of integrating over the probability simplex, whereby 

g(q)dq= n (^i-^jtW^'^'^^k =iVe-V^'^-^P(p)dpdr. (C.4) 

l<«<i<M k=l 
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The new variables qi take on values independently in the range [0, cxo). Integrating 
over all the values of the new variables, we find that the normalization constant is 
given by A^ = Q/T{^v), where Q = J Q(q)(iq. Similarly, we find that 



^/g^Q(q)dq = Q 



r(/i^) 



VMjPip) dp • 



(C.5) 



Notice that the first product in Eq. (C.4) is the square of the Van der Monde 
determinant [72, 73] 



A(q) = n (^^ - ^i) 

i<»<i<A' 



1 






1 
9a* 






rf{qi] 






l(9i) ••• ^m-i(9m) 



. (C.6) 



The second determinant in Eq. (C.6) follows from the basic property of invariance 
after adding a multiple of one row to another, with a = u — fi and the polyno- 
mials r^{q) = k\L1{q) judiciously chosen to be rescaled Laguerre polynomials [oti], 
satisfying the recursion relation 



iU 






and having the orthogonality property 



/ t/ge-Vr^(g)rng) = T{k + l)r(« + k + 1)5^ 
Jo 



(C. 
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These facts in hand, we can evaluate 



Q 









k=l ■ 



fi\Y[r{k)r{a + k) . 



(C.9) 



k=l 



For fi = u, a = and we can siniphfy the algebra considerably. Under these 
conditions, 



i,j=l -^ k,l=0 - 






Kl/2) Hl/2) _ / .(1/2) 



(CIO) 



where 



r(l/2) 



We thus have 



e '^y/qLk{q)Li{q) dq. 



1 ^"' 



fcJ=0 



.(1/2) .(1/2) _ /.(1/2) 



(C.ll) 



(C.12) 



except that the integral needs to be evaluated. For that, we use the generating 
function for Laguerre polynomials^ [36] 



;i-;z)-ie^"/^-^ = ^Li(x)2' 

1=0 



\z\ < 1, 



(C.13) 



and 



er-VL^it) it = ^''^ + I'*;'" 7; + 'K -^-'F (~„,3 + l;a + 1. 1 ) . (CM) 

nl L [a + 1) \ s 



^I thank C. Chandler for reminding me this trick. 
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F being the hypergeometric function such that 



oo 

F{a,b;c;z) = ^ 

n=0 



a)„ {b)n z"' 



c„ n! 



(C.15) 



and {a)n = a(a + l)(a + 2). ..(a + n — 1) is the Pochhammer symbol. Note that if a 
is a negative integer, (a)n = for n > |a| and the hypergeometric series terminates. 
Then, 



E 

1=0 






e-^v^Lfc(a;)(l - zy^e^^l^'^ dx 



s I e """^/xLkix) dx 



s = l/{l-z) 



r(3/2)r(A; + l) 

A;!rri) 



""'^•W) 



I - ^)*+V2 



^ Y- (-fe)t (3/2), 1 

^ (l)t ^! 
t=o ^ ' 

n^' i-iy{-k),{3/2), 1 



1=0 t=o 



(t!) 



(C.16) 



whereby 



r(l/2) 



^(-l) 



2 /! 



E 



:-k)ti3/2) 



(t!) 



-(t + b^ 



(C.17) 



and {a)n = a{a — l)(a — 2). ..(a — n + 1) is the 'falhng factorial'. Using the following 
identities for the Pochhammer symbols 



(x)„ = (-l)"(-x)„, 

-X)n = {-iy\x-n+l)n, 

{x)n = r(x + n)/r(a;), 



(C.18a) 
(C.18b) 
(C.18c) 
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we have 



;-i)' A / k \ [r(t + 3/2)]^ 



.(1/2) ^ l-^J V^ 

*=' /! ^ \ ^ / t!r(t-/ + 3/2) 



n i-iy ^ /3 3 



3-f^2 



4r[|-^]Z! " " V2'2' ' '2 



A;;l,--/;1 . (C.19) 



0.80 r 



0.78 



0.76- 



0.74 - 



0.72 



0.70 



10 



Figure CI: The value of the norinahzed negativity (A/')/2"'/^ of a random (Haar- 
distributed) pure state 



To get the final expression for the negativity in Eq C.12, we substitute the expres- 
sion for the integrals from Eq C.19. The expressions are not very illuminating, and 
for the lack of an asymptotic expression, we plot the numerical value in Fig C.l. As 
we see, the negativity also scales exponentially with the system size in the asymptotic 
limit. In fact, it goes as a constant multiple (0.720507) of the maximum possible 
negativity, just as the Schmidt rank of the DQCl circuit that has a random unitary 
in it. 
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